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ABSTRACT 


The successful prediction of the performance of a new or modified air- 
craft depends heavily on an accurate estimation of its lift and drag. The 
present work represents an effort to reduce the cruise drag of light aircraft 
through an analytical study of the contributions to the drag arising from the 
shape of the engine cowl and the forward fuselage area and also that resulting 
from the cooling air mass flow through intake and exhaust sites on the nacelle. 
It contains descriptions of the methods employed for the calculation of the 
potential flow about an arbitrary three-dimensional body with modifications 
to include the effects of boundary layer displacement thickness, a nonuniform 
onset flow field (such as that due to a rotating propeller) , and the presence 
of air intakes and exhausts. It also contains a simple, reliable, largely 
automated scheme to better define or change the shape of a body. 

A technique has been developed which can yield physically-acceptable skin- 
friction and pressure drag coefficients for isolated light aircraft bodies. 

For test cases on a blunt-nose Cessna 182 fuselage, the technique predicted 
drag reductions as much as 28.5% by body recontouring and proper placements 
and sizing of the cooling air Intakes and exhausts. 
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INTRODUCTION 


The successful prediction of the performance of a new or modified air- 
craft depends heavily on an accurate estimation of its lift and drag. Although 
the importance of these data is recognized , many light, aircraft manufacturers 
and most universities continue to depend essentially on the semi-empirical 
correlations of wind tunnel and flight test data along with personal rules of 
thumb in order to develop these estimates. 

For some time sophisticated techniques utilizing high-speed digital com- 
puters have been available to predict the performance characteristics of light 
aircraft. Although quite accurate, most of these techniques consume large 
amounts of computational time and computer core storage. Several attempts - 
some successful - have been made within the last five to ten years to reduce 
the time required to perform these estimates without significant sacrifices 
in accuracy. These technological. advances provide suitable starting points 
for the designer to "test" new or modified aircraft shapes without the expen- 
sive and time-consuming wind tunnel and flight tests. 

The present work represents an effort to reduce the cruise drag of light 
aircraft through an analytical study of the contributions to the drag arising 
from the engine cowl shape and the forward fuselage area and as well as that 
due to the cooling air mass. flowing through intake and exhaust sites on the 
nacelle. Since efficient fuel use is an increasingly important factor in 
general aviation operations, any design procedure which can lower the fuel 
consumption of a variety of aircraft through a drag clean-up is a welcomed 
advance in technology. 
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This report contains descriptions of. the methods. employed for the calcu- 
lation of the potential flow about an arbitrary three-dimensional body with 
appropriate modifications to include the effects of (a) the boundary layer 
over the body, (b) a nonuniform onset flow field about the body - typically 
produced by a rotating propeller, and (c) the presence of air intake and 
exhaust sites on the body for engine cooling purposes. The basic potential 
flow solution is accomplished by a solution of a Fredholm equation of the 
second kind, while the effects of the boundary layer, the nonuniform flow 
field, and the air intake and exhaust sites are included as the boundary con- 
ditions to the Fredholm equation. 

As an aid in the preparation of the input data to potential flow calcu- 
lations, this report also discusses a simple, reliable, largely automated 
geometry scheme to augment and/or to modify the body shape information by 
various techniques. Being expenditious and inexpensive, this scheme - a 
digital computer program - is a valuable tool to the researcher who wishes 
to better define or change the shape of a complete body or of regions on the 
body of particular interest before beginning the potential flow calculations. 
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THEORETICAL APPROACH OF HESS. AND SMITH FOR 
CALCULATION OF POTENTIAL FLOW ABOUT THREE- 
DIMENSIONAL NON-LIFTINO BODIES BY A SURFACE 
SOURCE-DENSITY DISTRIBUTION METHOD 


The problem under consideration is that of the potential flow of an 
incompressible, inviscid fluid about an arbitrary three-dimensional body. 

If the fluid density is constant and the viscosity is zero, the general 
Navier-Stokes equations reduce to the Eulerian equations of motion 

|^ + ( V • grad)$ * - ~ grad p ^ 

dt P 

where V is the fluid velocity at any point, p is the constant fluid density, 
and p is the fluid pressure. The continuity equation becomes 

div($) - 0 (2 


All body forces are assumed to be conservative and their potentials absorbed 
in the pressure Therefore equations (1) and (2) are valid expressions for 

the flow field exterior to the boundary surfaces. 

In order to discuss the flow about a three-dimensional body surface, let 
R* denote the exterior flow field and S denote the surface of the body (also 
the boundary of the region R). The body is assumed to have a surface repre- 
sented by an equation of the form 


S(x,y,z) - 0 


(3) 


where x, y, and z are the Cartesian coordinates as depicted in Figure 1. Under 
the assumption that the location of all boundary surfaces are known and that 



Figure 1. Body surface represented by an equation 
of the form S(x,y,z) ■ 0 


the normal component of fluid velocity is prescribed on these boundaries, 
the boundary conditions may be given as 

$ • n| a - F (4) 

where n is the unit outward normal vector at a point on S and F is a known 

function of position or time or both. To be complete, a regularity condition 

at infinity must be Imposed for the exterior flow problem. The onset flow $ 

00 

is assumed to be a uniform stream of unit magnitude; however, this restriction 
is not essential to the general derivation. 

it should be noted that the above equations do not define a potential 
flow. Since potential flow is a consequence of the condition of irrotation- 
allty* the usual approach to determine the equations of potential flow is to 
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assume that the velocity field $ is irrotational and therefore the velocity, 
field can be expressed as the negative gradient of a scalar potential function 

Letting the velocity field $ be the sum of the uniform, incompressible 
onset flow ^ and the perturbation velocity field v due to the surface 
boundaries, then 

V « V + v 


where 


v ■ - grad 4> 


( 6 ) 


Since the onset flow and the perturbation flow are jLncompressible , the con 
tinuity equation (2) is satisfied: 


div(V w ) - 0 
div(v) ■ 0 


div($) - 0 


As expected, the potential $ satisfies Laplace s equation 

V 2 $ • 0 


(7) 


in the region R’ exterior to surface S. By equation (4) the boundary condl 
tions on <|> become 


grad ♦ • n| 





F 


and the regularity condition for the exterior problem becomes 


(8) 


I grad <p | -► 0 at infinity 


C9 


Therefore equations (7), (8), end (9) represent the necessary potential flow 
equations to be solved . 

As demonstrated above, potential flow Is derived from the fact that ^he 
velocity field is determined by the continuity equation (2) and the condition 
of irrotationality (6). Thus equation (1) is not used, and the velocity may 
be determined independently of the pressure. Even though Laplace's equation 
(7) is the simplest and best known of all partial differential equations, the 
number of useful exact analytical solutions is quite small because of the dif- 
ficulty in satisfying the boundary conditions. Therefore indirect* methods of 
solution must be used to give satisfactory results from the various prescribed 
body surfaces and boundary conditions. 

The reduction of the problem to an integral equation for a source-density 
distribution on the body surface can be accomplished by the use of Green's 
theorem. The problem now is to reduce the potential flow equations (7), (8), 
and (9) to an integral equation. For a single three-dimensional body, con- 
sider a unit point source located at a point q whose Cartesian coordinates 
are x q » y^, and z^ in Figure 2. At a point P with coordinates x, y, and z, 
the potential due to this source is 


♦ 


1 

r(P,q) 


( 10 ) 


where r(P,q) is the distance between points P and q. The solution is 


Indirect or exact numerical methods contain the exact analytical formulation 
to the problem and have the property that the errors in the calculated results 
can be made as small as desired by refining the numerical procedures. Approx- 
imate methods contain analytical approximations in the formulation itself and 
thus places an accuracy limit on the results regardless of the numerical pro- 
cedures used. 
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Figure 2. Notation used in describing the surface 
source-distribution potential 


constructed of elementary potentials of the above form of an ensemble of 
sources. The resulting potential satisfies equation (9) and satisfies equa- 
tion (7) at all points except q. Because of the linearity cf the problem, 
such a potential that is internal or upon the boundary surface S satisfies 
equation (9) and satisfies equation (7) in the region R’ that is external to 
S. It is of considerable importance to determine the potential of a contin- 
uous source distribution on the surface S. 

Let o(q) be the local source distribution intensity, where the point q 
now represents a general point on the surface S. The potential of this dis- 


tribution is 
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By the procedure given by. Kellogg [3], the perturbation potential is given by 
equation (11) is differentiated, and the boundary condition (8) applied to it 
by permitting the point P to approach appoint. q on the surface S. The result 
is the following integral equation for the source-density distribution o(p): 


2tto( P ) * ^ In (r<MT) ° < q)dS * ‘ " ( P ) # K + F 


( 12 ) 


3 

where — denotes differentiation in the direction of the outward normal to the 
an 

surface S at the point p, and n(p) is the unit outward normal vector (written 
explicitly to show its dependence on location) . Equation (12) is a Fredholm 
integral equation of the second kind over the boundary surface S. 

The method of solution of Equation (12) is demanded to be numerical rather 
than analytical by the fact that the domain of integration is completely arbi- 
trary. The solution can be accomplished by first representing the body 
surface by a large number of small quadrilateral elements or "panels" (Figure 
3) . On each quadrilateral a control point is selected (usually the centroid) 
where the boundary condition is to be satisfied and where surface velocities 
are eventually calculated. A "matrix of influence coefficients", consisting 
of the complete set of velocities induced by the panels at each other's con- 
trol points, is then determined. The integral equation (12) is now approxi- 
mated by a set of linear algebraic equations for the values of the source 
strengths on the panels. Since each panel is assumed to have an independent 
value of constant source strength, the number of unknown parameters (source 
strengths) equals the number of panels or, more specifically, the number of 
linear equations. Once the source strengths are determined, the desired flow 
parameters may be calculated. This implementation renders the method as 
numerically "exact" and applicable to any arbitrary non-lifting body. 
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lOUNOmr LAYER SIMULATION FOR BODIES 
IN AN INVISCIO POTENTIAL FLOW FIELO 


Because of the considerations of large computation times, an ''exact 1 ' 
numerical solution of the complete Navier-Stokes equations is not feasible at 
this time. Yet the problem of determining the characteristics of a boundary 
layer flowing over a general, three-dimensional body is of great interest if 
real viscous flow effects are to be approximated. 

For essentially unseparated flow about a body at practical Reynolds 
numbers, the effects of viscosity are important only in a very thin boundary 
layer adjacent to the body's surface and in a thin wake downstream of the 
body. The following discussion will be restricted to the low-speed regime 
where the flow is incompressible or where the compressibility effects are 

sufficiently small to be handled by simple correction to an incompressible 
flow method. 

Smetana et al. [4] and Haaa [5] utilized a method of a two-dimensional 
boundary layer simulation technique expounded by Ll,hthlll [6]. This method, 
that of surface displacement or flow reduction, essentially consists of pro- 
ducing a thicker body by addin, the boundary layer dl.pl.ce.ent thickness to 
the original body in the direction along the local normal to the body surface. 
The potential flow about this modified body is thus the desired potentl.l flow 
to approximate viscous flow effects. This method is reasonably accurate if 
the body's cross-sectional ara. and volume do not change rapidly in the stream- 
wise direction, if no significant pressure gradients exist in the cross flow 
dlrsctlon, and if the wake i» adequetely modeled (see Sa.t.na £t al. [4], pp. 
125 - 132 ). 


10 


Since the present work Is essentially en extension of the work of Smetana 
et al, [4] of approximating real flows about light aircraft fuselages with an 
lnviscid, incompressible flow method, their procedure for the calculation of the 
displacement thickness, wall shear, and, ultimately, the lift and drag on a 
fuselage will be discussed briefly. It should be pointed out that a fuselage 
is generally considered to be a body with a plane or symmetry (normally the 
x-z plane) rather than an axis of symmetry. Effectively, this consideration 
simplifies the implementation of the method but also restricts its generality. 

Under certain conditions, as mentioned above, the three-dimensional 
boundary layer equations can be reduced to "simpler 1 * two-dimensional equations 
written in a general curvilinear coordinate system which describe the fuselage 
surface. By representing the fuselage locally by a section of a prolate 
spheroid with its major axis aligned with the local streamline* and its center 
in the x-z plane, the direction and magnitude of the flow velocity at one point 
on each panel may be determined by the induced velocity equations (6) . At the 
points where the boundary condition (8) is to be satisfied, these vectors (the 
streamlines) whose directions are the flow directions and whose lengths are 
proportional to the flow magnitudes can be easily determined. The method of 
isoclines may be used to sketch these streamlines; however, for computational 
purposes it is desirable to describe these lines in a more analytical fashion 
by assuming that for a given panel the flow direction across this panel is 
constant and that the quantity of flow is dependent on the distance and average 


*To adjacent streamlines form the boundaries of the flow of a given quan- 
tity of fluid. From the magnitude and direction of the flow over the surface, 
the position of these streamlines on the body may be determined. 
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velocity between adjacent streamlines. Therefore It la notad that alnca the 
quantity of fluid between two streamlines la always constant, the. flow velocity 
is Increased for converging streamlines and is decreased for diverging stream- 
lines. Also, this assumption is seen to be true only for Infinitesimal panel 
sizes. 

After fitting a section of the prolate spheriod to a section of the fuse- 
lage, as described by Smetana et al. [4] , the boundary layer equations written 
in general curvilinear coordinates are used with the local values of the sec- 
tion to determine the effects of body curvature on the displacement thickness 
and wall shear using a momentum integral formulation. At this point, the 
necessary means Jto describe the local streamwlse and crossflow coordinates on 
the body in terms of the reference coordinate system are available (as well 
as to write the boundary layer equations in a general curvilinear coordinate 
system) . 

To preserve the metrics of the general curvilinear coordinate system, the 
displacement thickness must be added normal to the body surface. In essence, 
a surface panel is translated parallel to its normal vector by an amount equal 
to the value of the displacement thickness calculated at the point where the 
boundary condition was satisfied. As shown in Figure 4, adjacent panels after 
the addition of the displacement thicknesses may not have coincident edges. 

This difficulty is quickly eliminated by simply averaging the new edge points 
1 and 2 to yield a new point. It is also observed in Figure 5 that a single 
line* dividing two strips of elements now becomes two lines. Again this dif- 
ficulty is remedied by averaging corresponding points on these lines to yield 

The connotation used above of "line" or "lines" actually represent a curve or 
curves in three-dimensional space. The sole purpose of this usage is the 
simplicity in viewing. 
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Figure 4. Displacement thickness addition to a body 

a single line. A word of caution should be mentioned at this point • It is 
recognized in Figures 4 and 5 that the averaging procedure can translate and 
rotate the panels. Therefore the metrics of the general curvilinear coordinate 
system are no longer preserved. The error introduced may be quite significant 
if the surface curvature varies rapidly from panel to adjacent panels. 

Since this method for determining the pressures and velocities over the 
surface of the body is an invlscid one, it always places a stagnation. point at 
the downstream end of the closed body. As a result, the method predicts stag- 
nation pressures at the aft end of the body in what is physically a wake region 
The pressures in this wake region is generally less than atmospheric. This 
relatively low pressure on the aft portion of the body as opposed to the high 
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Figure 5. Illustration of separated strips of elements 


pressure region near the nose resolves into a force acting in the flow direc- 
tion which is commonly known as form or pressure drag. It is apparent that 
for bodies producing regions of flow separation a reasonably accurate model 
to represent the wake effects must be used if meaningful drag results are to 
be obtained. 

Smetana et al. [4] replaced the physical wake by a solid extension of the 
physical body since the wake may be considered ss a region of "dead" air 
relative to the remaining flow field. They also assumed, rather arbitrarily, 







th.. this solid «t«n.ion (.» «k.-body) bo, In. .« .bout th. Uot t.n ».«.« 
of th. Ln.th.of th. phy.lc.l body, .nd «h.« «h. «.h.--b.dy .nrf.o. ,r.du.U, 

fair 8 via an expontantial function* 

Ih. .»rf.c. of thi. w.k.-body i. ...l«n.d qu.dril.t.r.1. ot p.«L of 

approximately th. .am. area *• tho.. on th. phy»io.l body, »• d.piot.d in 
Figure 6. 

OdXfIWt IOOY 


MCE IOOY 


Figure 6. Definition of wake body 

The prasaur.., a. d.t.mln.d by th. invl.cid flow method, on. tho.. p».l. of 
th. wake-body which 11. imm.dLt.1, .bov. th. .,uiv.l«t Ni on th. phy.lcal 
body .r. .ppll.d to th. p.n.l. on th. phy.lcl bod, .long th. norm.1. to th. 
physical body. Th. pr...ur.. ar. lnt.,r.t.d ov.r th. bod, .urfac, and th. 
reaultlng fore acting on th. phy.ic.1 bod, can b. n»d to find a lift and 
a drag. Sine. th. pr...ur«. on th. up.tr.am portion of th. wak. body will 
generally b. 1... than th... on th. ...r of th. phy.lc.l body ..cording to th. 
lnvlsdd flow c.mput.tlon, th. lnt.gr.tlon of for... will indie... . ».t dr.g 
on th. body. Th. tot.l dr.g on th. b.dy 1., of tour... th. .u. of th. pr...ur. 
drag and the akin friction drag. 
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Thi nathcd discussed above ia thataby characterised by the following atapa* 

(1) Tha aurfaca of tha leolated fuaalaga ia rapraaented by a sufficiently 
large number of quadrilaterale or four-eidad panala. 

(2) All four cornara of tha panel are moved into the aama plana through 
a procedure which determines the direction of the normal. 

(3) A source of undetermined strength is placed on each panel, and the 
prescribed normal boundary condition is required to be satisfied. 

(4) The resulting system of equations are solved for the source strengths. 

(5) The velocity over the body surface is calculated, ftnd^jth ^ str eaal ine s 
and surface pressures are determined. 

(6) Two-dimensional, momentum- integral-type boundary layer computations 
are performed along streamlines to find the local values of displace- 
ment thickness and wall shear. 

(7) The wall shear is integrated over the surface to find the skin fric- 
tion drag of the isolated fuselage. 

(8) The body shape is modified by attaching a wake-body toward the trail- 
ing edge and by accounting for the displacement thickness effects . 

(9) A new set of source strengths and surface pressures corresponding 
to the wake-body shape is calculated. 

(10) The surface pressures are Integrated to find the lift and pressure 

drag. 

(11) The total drag is determined from the sum of the skin friction drag 
and the pressure drag. 

This method does have imposed restrictions and limitations as discussed 
in [4]. It should be observed that the boundary-layer computations were allowed 
to iterate only once because of the great amount of work Involved in successively 


modifying the -body shop* to account for displacement thicknaaa effects. Tha 
displacement thickness ovar tha forward part of tha body is usually small and 
the use made of this information over the aft portion of the body in the present 
method's drag calculations is rather approximate. If the fuselage geometry is 
adequately presented to the inviscid flow field computation or if pressure 
gradients are such that flow separation is imminent only over the far aft por- 
tion of the body* the calculated results are quite acceptable and compare very 
well to experimental results. Otherwise* the user must examine the results 
with caution and, based on experiments and engineering judgement, decide their 
reliability. 
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PROPELLER NRKE SIMULATION 


Flo. vl.uallz.clon. of . typlcol nk. (ot .lip.tr.*> ■«*««• by « 
alrcr.ft prop.ll.r v.rlfy that prop* modeling of .uoh . Pben-n™ i. 
uthn.tlc.lly ntculy difficult, if not i«po..ibl.. for d.o.d.. InvMtl- 
gator. have ba.n .tudying thi. probln la the attnpt to develop adaouata 
uthnatical modal. to predict th. flo. ch.ract.rl.tlc. .ithin th. nk. *d 
the wake ' a influence on it. .unrounding.. Significant advancements haw b..n 
ud. through th. research on h.licopt.r rotor., .unuri..d by St.pni.nkl [7). 
Effectively, the *k. of h.licopt.r rotor, (in th. hovering nod.) uy be con- 
sldered to b. .ynonynou. .1th that of aircraft prop.ll.re, ...ning of cour.e 

that their differences are noted and understood. 

It 1. not th. intent of thi. report to expound on th. varlou. method. 

devised through th. decade, or to present a new ona. It 1. of concern to 
simulate the effect, of a propeller-lik. nk. upon the invi.cid Potential flo. 
calculation, of a thr..-di«n.ion.l body, unely, an aircraft fuselage or 
nacelle. Th. method to be utilized i. that cononly knon a. * pre.cribed, 

rigid wake.* 

in th. study of nk. .tructur.., vortar thsory 1. especially useful since 
th. observed b.h.vior of phy.lcl nk.. c.n b. msth.metic.lly .xpUlned through 
it. fundnsnt.l conc.pt.. Clanlcal exprn.ion. of th. Slot-S.v.rt 1.. can b. 
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applied whan astabllahlng tha relationship* batwaan tha atrangth and geometry 
of tha vortices, and tha velocities Induced by them In tha surrounding £1 » • 

Tha basic definitions and laws of vortex motion will be presently stated for 
later convenience. 

Even though one of the principal characteristics of ideal-fluid motion 
is its irrotationality, it is often necessary to consider flows in which a 
few points, lines, or regions of locally rotational flow. exit. It is often 
possible to analyze viscous fluid motions as being irrotational, except at 
the locations where the rotation is concentrated. 

The vorticity vector is simply defined to be twice the fluid rotation 
vector or to be the curl of the velocity vector. A vortex line is a line in 
the fluid giving the direction of the vorticity at every point on that line. 

A vortex tube is a tube of finite cross-sectional area on whose surface (known 
as a vortex surface ) the vorticity vector is always tangential. A vortex fila- 
ment is a "vortex tube" with an infinitesimal cross-sectional area, whose axis 
is a vortex line, and a finite value of circulation. 

The behavior of vortex filaments in an ideal fluid is governed by the 
following, theorems or laws of Helmholtz and Kelvin: 

1. The strength of vortex filament or tube is invariant at all cross 
sections along the axis. 

2. A vortex filament or tube cannot end in the fluid; that is, it muni, 
extend to the boundaries of the motion or form a closed loop. 

3. Tf an inviscid fluid, subject to conservative external forces, Is 
originally Irrotational, it will remain irrotational. 

As the tool for the determination of the induced flow field of a three- 
dimensional vortex system, the Biot-Savart law can be developed by using tin 
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concept of e vortex filament • Coneider en arbitrary filament L curving 
through space (Figure 7). The Incremental. velocity dv Induced at a point P 

f 



Figure 7. Velocity induction by element of 
three-dimensicuial ring vortex 


by an element dt of the filament L of circulation strength r is 


dv » (r/4ir)[(dt x ct)/d 3 ] . 
The. total .induced velocity at point P is 


<n) 


v - (r/4it)^ t (dt x t)/d 3 ] 


(l<4) 
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where the line integral indicates that the intagration is performed along the 
lina ..of the filament I* in the direction of dt with-a clockwise circulation 
around L, and 3 is the distance. .between the point P and the element dt* As 
indicated in the above figure and by the cross product dl x the velocity 
increment (13) due to the element is perpendicular to the plane of 3 and d£. 

Baskin et al^. [8] suggested the following procedure: The induced velocity 

vector (14) is resolved into the Cartesian components u, v, and w along the x, 
y, and z axes, respectively. Let the equation of the line L be given in the 
parametric form 

5 « 5(8) , n ■ n(0) » 5 ■ C(0) 


where 8_is the .parameter (being suggestive as an angle for the curved filament). 
As the parameter 6 varies from its initial value 8^ to its final value 8^* 
point Q(5,n»0 describes the curve L. The vectors 3 and d t can be expressed as 


3 - (5 - x)i + (n - y)j + (? - z)k 

d£ - ( H de)i + de)3 + (^ de)k 


(16) 


where i, j. and k are unit vectors of the x, y, z coordinate system. 

Substituting the above expression into (14), the x, y, and z components 
of the induced velocity are determined to be 


» • fc e [ £ - «> - ■ ">$ (1,) 

v ‘fe I l39 (x ' {) * § <! ' {,1 7 <18) 
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where 


* • fc a | £ ld8 <y * "> * d? tx ' S) 1^3 

d * /(x - O 2 + (y “ n) 2 + ■(* * ?) 2 • 


( 19 ) 


( 20 ) 


Emanating from the blades of the propeller, free vortices form the wake. 
These free vortices can be subdivided into shed vortices which, at the moment 
of leaving the blade, are. parallel to its axis, and trailing vortices outflow- 
ing along the blade span in the direction either perpendicular or. approximately 
perpendicular to the blade axis as depicted by Figure 8. Among the trailing 
vortices, the tip vortices (those leaving the blade tips) usually dominate the 



Figure 8. Formation of the wake by free vortices 
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flow picture in all regions of operation. In the complete vortex system the 
shed and trailing vortices may form a surface of vorticity of a generally 
helical shape behind each blade as indicated by Figure 9. 



Figure 9. Helical-shaped surface of vorticity 


It is easy to comprehend that, in reality, the geometry of the wake and the 
strengths of the vortices forming that wake may vary with time because of the 
interaction between the vortices, flow fluctuations, blade rotation, etc. 

In this report a simple concept* of the wake will be used since the sole 
purpose of the simulation is to superimpose the resulting nonuniform induced 
flow field on a three-dimensional body for an inviscid potential flow calcu- 
lation. Since the tip vortices are most prevalent, it seems natural to assume 


*More sophisticated physicomathematical methods for modeling the wake structure 
are available for computer implementation. These techniques, as applied to 
the analysis of helicopter rotors, are discussed by Stepniewaki [/]. 
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that the circulation around the blade remains nearly constant along its entire 
span and to approximate the wake structure by. modeling only the tip vortices*™. 
The tip vortices will trace out roughly the paths traveled by the tips of the 
propeller blades. For a two-bladed propeller. Figure 10 represents the system 
of trailing vortices of a propeller with a constant circulation from the root 
to the tip of the blades. 


r 



Figure 10. System of trailing vortices of a propeller 
with a constant circulation 

These vortex lines or filaments constitute the slipstream (wake), and the fluid 
motion in this slipstream can be determined as the induced velocity of the 
system of vortices. The fluid in the slipstream has an increased axial velocity, 
and the rotational velocity is in the same sense as the rotation of the propel- 
ler. Therefore, the general behavior of the flow is the same as that postulated 
by the momentum theory. 

Physically, the strength of the circulation will vary along the blade. 

Due to this variation, trailing vortices will arise, not only at the root and 
tip of the blade, but from every point of its trailing edge as shown in Figures 
8 and 9. Since the induced velocity of a system of helical-shaped vortex sur- 
faces which constitute the propeller's slipstream is difficult to calculate. 
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the usual procedure Is to assume that the propeller has a very large number of 
blades or that the propeller -has a near-zero advance angle. This assumption 
implies that the vorticity of the slipstream is distributed throughout the 
system instead of being concentrated on a small number of vortex surfaces. 


Instead of the helical surfaces, the slipstream may be considered as a 
close succession of vortex rings. These vortex rings can be applied to and 
within the boundaries of the slipstream as shown in Figure 11. 



Figure 11. Application of vortex rings 


The vortex rings applied at the boundaries of the slipstream represent the tip 
vortices while those applied within the boundaries represent the trailing 
vortices springing from the blade due to the variation of circulation along 
the blade. The vorticity of such systems of vortex rings represents the in- 
creased axial velocity of the slipstream, and the vorticity of the rectilinear 
vortex, constituting the sum of the root trailing vortices, represents the 
rotation of the slipstream. The whole slipstream can be considered to be full 
of vortex systems of this simple type. 

If the slipstream is allowed to contract as it would in real flows, the 
condition of continuity of the flow requires that the axial velocity and the 
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angular velocity must increase as the radius of the slipstream decreases. If 
the contraction of the slipstream can be ignored, the axial, velocity in the 
ultimate wake must be equal to twice that at the plane of the propeller.* 

This result is exactly that found by the momentum theory for propellers [9]. 

According to the theorems of Helmholtz and Kelvin, there is no exchange 
of either mass or momentum between the vortex filament and the rest of the 
fluid. Hence, if a vortex filament were located in a mass of moving fluid, 
it would move with the fluid. Thus velocity fields induced by a system of 
vortices can, in t urn . produce motion of those vortices belonging to the 
system. 

At this point it is informative to discuss the motion of ring vortices 
if the vortices are free to convect in the fluid. For a single ring vortex 
(Figure 12a), the streamlines of the induced velocity (Figure 12b) may be 
divided into their velocity components and vectorially summed to obtain only 
an axially-directed velocity, since the y- and z-components sum to zero. 
Therefore, a single ring vortex translates along its axis due to its own 
induced velocity field. 

The motion of two ring vortices of equal strength, size, and sense on the 
same axis is easily comprehended by considering their mutual influences. Each 
ring induces components of velocity into the core of the other as sketched for 
a meridional plane in Figure 13. The forward ring receives an outward induced 


*With no slipstream contraction, the vortex system is simply a long cylinder 
of vorticity extending indefinitely in one direction from the propeller disc. 
The induced axial velocity at a point in the wake, where the cylinder extends 
Indefinitely in both directions, must be double that at the corresponding 
point in the propeller disc. 
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Figure 12. Motion of a single ring vortex 
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Figure 13. Mutual influences of two ring vortices 


velocity component, while the rearward ring receives an inward component. As 
its diameter increases, the forward ring decreases in speed, according to 
equation (14), while for the same reason the rearward ring increases in speed 







srsss#* 

due to its decreasing diameter. Soon the smaller , faster rearward ring passes 
through the larger, .slower forward ring x „The roles_of the vortices are now 
reversed. The smaller ring enlarges and its speed decreases, while the larger 
one decreases in diameter and speeds up. This leapfrog process, sketched in 
Figure 14, continues endlessly in an inviscid fluid. 


Figure 14. Leapfrog action of two free vortex rings 


In the present method the propeller is assumed to have a large number of 
blades. The approaching flow to the propeller is to be steady and uniform, 
while the receding flow within the propeller wake is to be nonuniform. The 
wake is also assumed to be a three-dimensionally rigid and prescribed one 
consisting of a succession of ring vortices representing only the trailing 



vortices, with the exception of the root trailing vortices, emanating from 
along the span of the blades. The root trailing vortices, which are simply 
rectilinear vortices coexiting along the propeller axis, would be the primary 
contributors to the angular motion of the slipstream, and there are basically 
two reasons for not including them in this method. The first and foremost 
reason is that the boundary layer method described earlier does not allow for 
crossflows on the body. Secondly, flow visualization studies indicate that 
the self-induced effects of the vorticity rapidly, carry the root vortices up 
through the center of the slipstream where they are dissipated. It is further 
assumed that the remaining trailing vortices quickly roll up within a few 
chord lengths from the plane of the propeller to form a single concentrated 
vortex trailing, from near the tip of each blade. This analytically-simplify- 
ing assumption is also verifiable from flow visualization studies. 

Since the slipstream is to be rigidly prescribed, certain aspects must be 
considered for the proper location and strengths of the vortex rings. As 
stated earlier, the trailing vortices quickly roll up to form a single trailing 
vortex. It therefore^ seems pausible that in the plane of the propeller several 
concentric ring vortices, each of different radius and strength, may be speci- 
fied, while fewer and fewer vortex rings may be specified further downstream. 
The location of the downstream vortices may be specified rather arbitrarily*, 
but their strengths should satisfy some mathematical requirement or some phys- 
ical phenomenon. The laws of Helmholtz and Kelvin state that the circulation 
is conserved. In other words, if a plane of the slipstream contains two vortex 


*An attempt should be made to approximate the results of momentum theory for 
propellers: i,.e. , the induced velocity in the ultimate wake be twice that in 
the plane of the propeller. . 
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rings of circulation strengths and r 2 while another plane further down- 
stream contains only one ring, then the strength of the single ring of the 
downstream plane must be + r 2 » Under the present method* the strengths of 
the vortex rings may be arbitrarily specified* with the differences between 
their theoretical strengths and those actually used being accounted for by 
assuming some type of viscous dissipation. 

The effect of a ring. vortex at a sufficient distance from itself is 
approximately the same as that of a three-dimensional doublet of strength 
1/2 Tr 2 , where r is the strength of the vortex ring and r is its radius [lO]. 
This statement results from the observation that as the distance from the ring 
vortex becomes greater and greater, the details of the ring's shape become 
less and less important. Since reductions in computational times are always 
desirable, the present method incorporates doublet flow whenever possible. 

The computational savings are quite substantial since the calculation of the 
velocity induced by a doublet does not require integration. 

In order to apply doublet flow, its theory will be discussed briefly: 
Consider the situation where a source of strength q and a sink of equal 
strength are located a very small distance i apart. Letting the distance l 
go to zero while q is allowed to increase indefinitely in such a way that the 
product q& remains finite and equal to a constant v* then the potential at a 
point P due to this so-called doublet at s is given by [Figure 15]. 
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DOUBLET 


Figure 15. Doublet potential definition 


where y is strength of the doublet in the direction of its axis. The positive 

doublet axis is simply the direction from the sink to the source. 

* 

The velocity field induced by a doublet is readily obtained from the 
equation 


$ - grad 4> - - ^ grad 


u*r 


The velocity components in Cartesian coordinates are thus 
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where m*. w . end 4 « g are the doublet’s strength in the x, y, * directions, 
respectively, and r x , r y , and r g are the x, y, * components, respectively of 
the vector r (Figure 15). 

The substitution of doublet flow for ring vortex flow requires that the 
origin of the doublet be the same as that of the ring vortex, the doublet axis 
be perpendicular to plane of the ring vortex and is directed so that the 
induced velocity is of the same sense as that of the ring vortex, and the 
doublet strength y be a "suitable” function of the ring vortex strength F. 

At sufficiently large distances from the orig’n of a doublet of strength 
1/2 Tr 2 , the induced velocity of the doublet approximates that of a ring vortex 
with its origin the same as that of the doublet. Unfortunately this doublet- 
strength is only valid at large distances. In the present method, the doublet 
strength 1/2 Tr 2 was reduced by a factor of 0.00932 so that the velocity in- 
duced by the., doublet closely matches that of a ring vortex of strength r at 
and beyond a distance of 1.5 times the diameter of the vortex ring (Figure 16). 
For all points at a distance greater than 1.5 times the vortex ring diameter, 
the induced velocity is calculated by the doublet flow equations. Otherwise, 
the vortex ring equations are used. 

Since vortex filaments are discontinuities in an inviscid fluid, it is 
important that no filament intersects with the body's surface. Otherwise, 
extremely large induced velocities exist on the surrounding surface, as pre- 
dicted by equations 13, 14, 17, 18, or 19 since d •* 0. In determining the 
potent .1 flow solution of a body in the slipstream of a propeller, the 
pres method calculates the diameters of the vortex rings so that no vortex 

in*-- ts with the body and that the mass flow between the body and the 
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DISTANCE FROM ORIGIN OF DOUBLET OR RINO VORTEX 


Figure 16. Matching doublet flow to vortex flow 

boundary of the slipstream is constant. These requirements are sufficiently 

satisfied by the following procedure: 

Each axial station of the body is described by a set of Cartesian points, 
each point havin 8 an identical x-coordlnate value (Figure 17). In other words 
an axial station is contained in a plane perpendicular to the x-axis. Assum- 
ing that each body cross-section is symmetrical about the x-z plane and con- 
tains the same number of points, only half of the total number of points must 
be supplied. The body must be closed at both ends, i.e., the first and last 
axial stations must not only contain the same number of points, each point of 
the same station must be equal. 

For each axial station, the cross-sectional area is computed by summing 
the individual areas of plane triangles of the cross-section. Each plane 
triangle is defined by the origin of the cross-section and two consecutive 
points on the body’s surface. The maximum distance and the average distance 
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Figure 17. Denotation of first and last body axial station 
with illustration of point numbering 


are also computed. The maximum distance is that from the origin to the 
farthest surface point, while the average distance is considered to be the 
radius of a circle equal in area as that of the cross-section of the body. 

With the maximum and average distances and the cross-sectional areas now known, 
an equivalent circular body is produced by letting the average distances be 
the radii of area-equivalent circles. This assumption allows for a crude, 
but sufficiently accurate*, method to satisfy the mass flow requirements and 
to guarantee that the vortex rings do not intersect with the body. 


The present method would not produce satisfactory flow fields for wide 
flat bodies with the equivalent circular-body concept. For these cases the use 
of an equivalent ellipsoidal-body concept along with elliptically-shaped vortices 
would be more suitable and realistic. 
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Consider an equivalent circular body within a slipstream (Figure 18a) and 
a typical cross-section (Figure 18b). The mass flow at the beginning of the 
slipstream Is given by 

prd 2 V/4 

where p is the density of the fluid, V is the mass-mean velocity of the fluid 
across the cross-section, and d is the diameter of the circular slipstream at 
its beginning. Similarly, the mass flow in an. annular ring is given by 

t 2 2 \xt 

pir(a - z )V 

where a is the distance from the body centerline to the boundary of the slip- 
stream, and z is the distance from the body centerline to the body surface. 
Requiring that these two_mass flows to be equal, it is found that the diameter 
of the vortex ring (or the boundary of the sli pstre am) must be given by 

a - A 2 /A + z 2 (21) 

As the body radius z varies, the diameter of the vortex ring varies accordingly. 
It should be noted that this is another way of saying that the radius and 
circulation were changed in such a way that the velocity induced on the body 
surface is the same as that at z * 0 i£ the body were not present. 

The origins of the various ring vortices may be assigned arbitrarily 
along the centerline of the body. Therefore, appropriate averages of the 
relative distances and areas must be calculated not only to suitably fair the 
given cross-sections, but to determine the diameters of the vortices. 
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Even though equation (21) eatisfies the continuity requirements at each 
position along the body, it does not guarantee that the vortex rings do not 
intersect with the body because of the assumption of the equivalent circular 
body. To eliminate this problem, a comparison is made between-ihe -diameters 
calculated by equation (21) and the maximum diameters of the original- body. 
The largest diameter is always used. At this point it is clear that if the 
cross-sections of the original body is greatly different from circular ones, 

the present method may not be sufficiently accurate for the mass flow calcu- 
lations. 

The method also has an automated procedure to generate additional vortex 
rings about and aft of a body by using exponential functions to assign the 
locations of their origins. The exponential functions are in turn expressions 
of the varying body dimensions. These functions may be terminologically 
defined as "weak forcing function" since their detailed behaviors are body- 
dependent, while their broad overall behaviors are set a prior . Basically, 
the method contains two such functions: one revelant from the start of the 

slipstream to the plane of maximum body width, and the other revelant from 
the maximum body-width plane to the assigned maximum axial extent of the 
vortices. Figure 19 depicts a vortex ring distribution and diameter variation 
about an aircraft fuselage in the wake of a propeller. 

As may be evident from the discussion above, accurate modeling of a 
propeller’s wake or slipstream is a very difficult task. Such a physico- 
raathematical model would necessarily be unduly elaborate (for the use intended 
here) and its computational requirements enormous. Since a major objective 
of the present work is to simulate the effects of the nohuniform flow field 
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Figure 19. Vortex-ring distribution and diameter variation 
about an aircraft fuselage in the wake of a propeller 


of a propeller-like wake within the context of the inviscid potential flow 
calculations over a three-dimensional body, accurate modeling does not seem 
to be warranted. Although no claims of accuracy are professed for this wake 
analysis, it would seem to provide at least some insights into the flow 


behavior within a slipstream. 
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FLOM PORT LOCATIONS FOR MINIMUM DRAG 
AND EFFECTIVE ENGINE COOLING 


Properly-located air Intake and exhaust sites provide adequate .cooling 
and carburetor induction and minimum pressure loss. Since the pressure loss 
(or "cooling” drag) is larger than necessary in many general aviations designs, 
the present method can be directed to determine cowl air intakes and exhaust 
sites that produce favorable interactions with the external flows and also 
provide good entrapment and extraction of cooling air. Favorable interactions 
with external flows may be interpreted as those where the inlet will capture 
sufficient cooling and carburetor air mass with a minimum disruption to the 
external flow and the exhaust will provide adequate extraction and contribute, 
if possible, to the streamlining of the aircraft. 

In this analysis the flow is assumed to be . incompressible and its free- 
stream velocity equal to unity. It is also assumed that a suitable schematic 
of the cowl interior of a body with inlet and exhaust flow may be given by 

Figure 20. In this figure and in the following equations,. Cp and Cp denote 

X XX 

interior pressure coefficients and EOA denotes the effective orifice area. 

An interior pressure coefficient may be defined as the pressure coeffi- 
cient required to balance the mass flows between the inlet and exhaust. In 
the absence of an engine, the interior pressure coefficient is constant every- 
where throughout the cowl interior. Since an engine occupies a large region 
of the interior, it obstructs the mass flow with its presence felt throughout 
the interior. A pressure drop will therefore arise between the pressure just 

inside the flow inlet and that just inside the flow exhaust as a result of 
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Figure 20. Schematic of cowl interior 

viscous losses incurred by the fluid as it moves through the engine fins and 
around baffles. With an engine present, two "different" interior pressure 
coefficients are thus necessary to include the effects of engine fin and 
baffle pressure drops. 

The effective orifice area depends directly on the engine geometry and 
the ''tightness'' of the engine inside the cowl. The tightness of an engine 
inside the cowl refers to the spacing between the cowl and engine. If this 
space exists, air mass will tend to leak around (or about) the engine instead 
of flowing through its fins and baffles and, therefore, reducing the cooling 
capacities of the air flow. With the assumption of a tight engine, the 
effective orifice area is that ''area” seen by the mass flow through the fins 
and baffles to provide engine cooling. Without the assumption, the effective 
orifice area must be some compromise of the leakage and cooling flows. 

From Figure 20, the mass flow m at each opening must be given by 
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where p, A, and V are the fluid density, the cross-sectional area, and the 
fluid velocity normal to the opening, respectively. Since mass is conserved 
the mass flow at each opening must be equal. That is, 


m * p 


'in 


in 


EOA 


EOA 


out 


out 


constant 


(23) 


By Bernoulli's equation for incompressible flow, 


». + | ■ ■> + 2 ^ 


(24) 


where p and V denote the reference freestream pressure and velocity, 
r 00 00 

respectively. Rearranging (24) and dividing by the reference dynamic 
pressure (1/2 pV 2 ) yields the pressure coefficient 


P " P. 

Cp - 2 

1/2 P< 


(25) 


For incompressible flows, the dynamic pressure is simply a difference of 
two pressures . Thus , 

Dynamic pressure - 1/2 pV 2 « Ap - P 2 - P x - pressure difference 
or 


V 


2 




(26) 


Subtracting p^ from Pg and adding p w to p^ of (26) yields 

V 2 • f [<P 2 - P.) - <Pi - p„>] 

2 

Multiplying and dividing the right-hand side of (27) by produces 
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or 



<p 2 “ 

1/2 P»* 


« 

1/2 pV* 




- Cp 


« 

1 


v- v. /Cp 2 - C ?1 


(28) 


Thus, equation (28) relates the velocity normal to an opening with the pressure 
coefficients on each side of the opening. 

At this point it is important to recapitulate on the representation of 
the body and to define parameters pertinent for the continuance of this method. 

As stated in the previous sections and shown in Figure 3, the three-dimensional 
body is constructed of an orderly arrangement of plane quadrilaterals or panels. 
This body representation facilitates the assignment of the inlet and exhaust 
sites because any panel or combination of panels may be declared as either 
inlet -openings or exhaust openings. Since stagnation flows are not considered, 
a body with an inlet site must, in turn, have an exhaust site. 

The total inlet area A. is assumed to be the sum of the individual inlet 

in 

panel areas, while the total exhaust area A o ^ is assumed to be the sum of the 
individual exhaust panel areas. Since the interior of the body (or cowl) is — 
treated essentially as a black box, it is necessary to propose a single inlet 
pressure coefficient Cp in and a single exhuast pressure coefficient CP 0Ut * 

These pressure coefficients are surmised to be area-averaged values of the 
revelant pressure coefficients of the selected inlet and exhaust panels from 
a solid-body (no inlet or exhaust) potential flow solution. That is, 


42 


ORIGINAL PAGE IS 
OF POOR QUALITY 


and 



Inlets 


1 A t Cp 

i=l L_ 

Inlets 
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Cp 


Exhausts 
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i-1 
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A l cp l 


Exhausts 
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The solid-body pressure coefficients are implied to remain equal for the 
opened body. This implication is only valid if the inlet and exhaust_panel 
geometries are relatively small compared to the overall body geometry. For 
this reason a "grid" refinement may be necessary to accurately model or shape 
the body particularly in the neighborhood of the inlet/exhaust areas. 

It is now possible to rewrite (23) by using (28). Hence, (V^ - 1) 


- CP. - pEOA/Cp, - Cp„ - PA„„ t /Cp„ - Cp 


XX out XX out 


or 


A lt /Cp ln - Cp. - EOA/Cp. - Cp. v - A„ llt /Cp„ - Cp 


XX out XX r out 


(29) 


To eliminate the square root, each "term" of (29) is squared; 


© 


© 


© 


4. <Cp in * Cp *> ' EOa2(Cp * - Cp **> ■ A ou t (Cp x* ' Cp o u t> 


(30) 


The rearrangement of terms and of equation (30) yields 
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4 Cp l» ' Cp x (A ln + ‘ * 


(31) 


while the rearrangement of terms ^2) and yields 


- A out Cp °ut + Cp xx (A out + B0 * 2) ' EOa2cP » 


(32) 


Upon solution of squstions (31) End (32) » 


CP 


xx 


A out Cp out + .2 


EM 2 • A 2 n • Cp ln 


A in + E0A 


A* + EOA 2 
In 


A L • A out + E0A2(A ln + A o»t> 


(33) 


and 


+ E ° a2<:P xx 

Pjt A 2 + EOA 2 

In 


(34) 


Equations (33) and (34) satisfy the limiting conditions of the effective 
orifice area: 


(1) 

EOA “ 0 

Cp - 

XX 

°»out ■ 

(2) 

EOA ® 

cp : 
r xx 

CP* • 



CP 


in 


The first condition implies complete flow blockage or no mass flow, while the 
second condition implies essentially unrestricted mass flow. 

Since the engine is far less than 100 percent efficient, it is necessary 
in a constant density representation to add to the cooling flow a quantity 
representing the internal heat addition to the air mass. This quantity is 
assumed to be physically represented by the difference between the maximum 
energy content of the fuel being used and the actual energy converted to 
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kinetic energy (developed engine power). The energy difference is surmised 
to be transferred to the surroundings dlrsctly in the form of heat which in 
turn leads to an Increase in temperature. 

The heat quantity Q is equal to the change in enthalpy H; i_.a. , 

Q ■ AH (35) 

Assuming that air is an ideal gas and the temperature variation is relatively 
small. 


AH ■ mCp air AT 


(constant-pressure process) 


(36) 


where m is the air mass, Cp a ^ r is the constant-pressure specific heat of air, 
and AT is the temperature difference. In terms of their rates, equation (35) 
becomes 


while equation (36) becomes 


Thus 


Q - AH , 


AH - mCp , AT . 
air 


• • a 


Q - AH - mCp aij .AT 


or 


AT 


. — 2 . 


mCp 


air 


(37) 


Dividing (37) by the reference temperature T yields 


T - T 

AT „ l _m 

T T 


“ Cp air T » 


(38) 
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The heat transfer $ ia now assumed to be 

$ ■ Hj • ffr - P 


where H f is the heat of combustion of the fuel, ffr Is the fuel flow rate, 
and P is the developed engine power. The fuel flow rate, though, is the ... 
product of the specific fuel consumption c and the developed engine powers 

ffr - c • P 

The mass flow rate m is by definition equal to equation (22) where p is taken 
to be the freestream density, A the total inlet area A in » and V the velocity 

normal to the inlet (- V/Cp in - Cp x ) . Therefore, equation (38) now becomes 


or 


T - T P[H, • c - 1] 

oo r 

T - P.ValrV^ - CP * 



P[H f • c - 1] 


) T Cp . v /cjT - Cp. 
oo oo *air °° in 3 


+ 1 


(39) 


Air has been assumed to be an ideal gas, and hence any increase in its 
temperature decreases its density. Since the incompressible theory on which 
this method is based assumes a constant density, the internal heat release 
must be represented by a compensating increase in the flow velocities in order 
to balance the mass flows. This compensation may be accomplished by varying 
the effective orifice area which, consequently, varies the effective flow 
velocities. The effective orifice area variation is determined as follows! 
The mass flow rate in the effective orifice region is given by (22) where A 
is the effective orifice area EOA and V is determined from (26) to be 
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( 40 ) 


The density change with temperature 
of state to be 


for an ideal gas is given by the equation 


P 1 T . 

— — m — 

p- T i 


or 


P 1 “ P * T. 


(41) 


A "new" mass flow rate 
(40) to be the density of (41): 


i is calculated by assuming the density of (22) and 


m. 


p. 5=] ™ V 


2hp 


<P. T .'V 


or 


“1 ■ { p « E0A Vif} V - 


(42) 


The mass flow rate (42) is seen to be the product of the otigi-k orifice flow 
cate and the square root of the temperature ratio T./ V Associating the 
temperature ratio with the original orifice area produce, a new effective 

orifice area 

eoa hew ■ E0A V T 
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Applying E 0A^g W to equations (33) and (34), new interior pressure coefficients 
are calculated. With T. > T (heat added), 

i OS 

eoa new < E0A * — 

Consequently, the flow velocities must vary in a constant density analysis in 
order to balance the mass flows. 

The typical behaviors of the Interior pressure coefficients due to vari - . 
ations in the freestream density, the inlet velocity, the freestream tempera'** 
ture, or the effective orifice area are shown in Figure 21a, while those due 
to variations in the fuel’s heat of combustion, the specific fuel consumption, 
or the developed engine power are shown in Figure 21b. 




Figure 21. Typical behaviors of the interior pressure coefficients 

with parameter variation 

On a percent basis, a variation in the effective orifice area or the specific 
fuel consumption had the most pronounced effect on Cp and Cp . Not only 
does Figure 21 indicate the slopes of the "lines" representing Cp and Cp , 

X XX 


it should be noticed that e measure of the "driving potential" of the air mass 
can be observed. This "driving potential" or pressure difference of Figure 
21a is seen to become smaller with increasing values of the pertinent param- 
eters, while that of Figure 21b is seen to become greater with increasing 
values of its pertinent parameters. Consequently, a measure of control over 
the quantity of air flow is possible by a strategic selection of the parameters. 

It is now possible to calculate the normal velocities on the specified 
inlet and exhaust panels and apply them as the complementary boundary condi- 
tions for the potential flow calculations. The normal velocities on the inlet 
and exhaust panels are surmised to be functions of the Interior pressure 
coefficients and the pressure coefficients that existed on those panels when 
the body was closed (with no internal mass flow). Thusly, 

v i ■ v . /c Pi - c »x (43> 

for the inlet panels and 

V - V /CP - Cp, (44) 

i * r xx i 

for the exhaust panels, where Cp.^ is panel i’s closed-body pressure coefficient 
and V (■ 1) is the freestream velocity. Attributable to the proposition of a 

OP 

single inlet pressure coefficient Cp in and a single exhaust pressure coeffi- 
cient Cp in the determination of Cp and Cp , the discriminants of (43) 

OtlC X An 

and (44) is. occasionally negative. This difficulty of imaginary numbers is 
eliminated by proposing that (43) and (44) are essentially 


v i - v » / 1cp7 - c p x i 


(45) 


T 


T 
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for the inlet panels and 

V t - ». -T^r^TT <46> 

for the exhaust pane!.. Th. normal velocities .« <«) - <«> «« — * 
to bo the complementary boundary conditions tor the potential flow solution 
discussed earlier. Th... normal velocities are Inserted Into (8) a. the 
now-known function F. 

in summary, the mas. flow of air through a body 1. dependent In this 
treatment upon (a) the pressures which existed when the Inlet end exhaust 
panels were closed, (b) any Internal pressure drops due to engine fin. and 
baffles, (e) the Internal heat addition, and (d) the are. ratio, of the Inlet, 
the effective orifice, and the exhaust. The present method consist, basically 

of the following steps: 

(1) Panels are chosen to represent the inlet and exhaust sites. 

(2) The total inlet and exhaust areas are determined by summing the 

individual areas. 

(3) An area-averaged pressure coefficient Is calcul.ted for th. inlet 
site and the sxhaust site.. 

(A) Internal heat addition is invoked through appropriate changes in 
the effective orifice area. 

(5) The interior pressure coefficients are determined, and the normal 
velocities on each individual inlet and exhaust panels are calcu- 


lated. 


A GRID REFINEMENT SCHEME FOR MODIFICATION 
OF A BODY'S GEOMETRIC DATA 


Since the preparation or modification of the geometric data of a three- 
dimensional body for input into a potential-flow program is tedious and time- 
consuming) a scheme has been developed to simplify this task. Given a set of 
Cartesian_-coordinate data descriptive of the body, this scheme may be used to 
correc t body misrepresentations , to change the body geometry , to refine the 
network (or grid) of the panels or .quadrilaterals that form the surface of the 
body, and to plot various orthographic, perspective, and stereoscopic views of 
the original and the modified body. 

The general procedure for specifying (or inputting) body points is the 
same as that given in Reference 4 for the NCSU BODY program. For clarity, 
this procedure is now reviewed: 

The body's surface is constructed by an arrangement of quadrilaterals 
with their corner points forming a network of intersecting lines to be called, 
from this time on, M-lines and N-lines (Figure 22). The M-lines are those 
connecting. corresponding points on the N-lines and are generally those running 
along in a direction near to that of the major body axis, while the N-lines 
are those surrounding the perimeter of the body in parallel planes perpendic- 
ular to the major body axis. Every N-line or M-line must have the same number 
of "defining" points as every other N-line or M-line with the stipulation that 
no N-line crosses another N-line and no M-line crosses another M-line, 
although they may converge to a common point. 

Since the primary purpose of the scheme is to generate data compatible to 
the NCSU BODY program and the present report’s potential-flow program which 
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Figure 22. Schematic of indexing scheme used for a 3-1 ellipsoid 
with 40 panels describing the half-body 
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assumes the body is symmetrical about the X-Z plane, only half of the body is 
necessary for input. For compatibility, the first N-line must be that at the 
front of the body (first to encounter freeetream air flow) with the remaining 
N-lines numbered consecutively from front to back. The points on each N-line 
are ordered such that the first point lies in the lower portion of the X-Z 
plane at the "bottom" of the body (Figure 22, Section A) with the remaining 
points consecutively increasing counter-clockwise (looking aft) to the upper 
portion of the X-Z plane at the "top" of *~he body. These points, i.e. , the 
corresponding points on successive N-lines, generate the M-lines. 

The workhorse of this scheme is a cubic-spline curve-fitting method [ll] 
coupled to a coordinate-system rotation-translation technique. Although cubic- 
spline fits are generally considered to be the smoothest of all curve fits 
[12], they often have difficulty in regions of extremely high curvature, where 
oscillations become magnified. Splining points against arc length alleviates 
the oscillations; the present investigators, however, chose to remove the 
oscillations by rotating and translating the body points in the reference 
coordinate system to points in the new coordinate system. 

As an illustration of the rotation and translation of the body points, 
consider the upper M-line* of the X-Z plane for some typical three-dimensional 
body (Figure 23). A curve fit of the points 1 through 23 - with their values 
given with respect to the reference coordinate system XYZ - may not be satis- 
factory between points 1, 2, and possibly 3 due to the presence of the high 


*A three-dimensional body contains two M-lines in the X-Z plane - an upper 
one and a lower one. Since M-lines are increasingly numbered consecutively from 
the body's bottom to its top, the upper M-line corresponds to the maximum- 
number M-line whereas the lower M-line corresponds to the minimum-number M-line. 
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slope between points 1 and 2. The alternative procedure is to transform the 
body points to those appropriate to the new coordinate systems. For instance, 
for the transformation of points 1 through 12 to those in the X’Y’Z* system, 
the following equations may be used:. 


where 


X i " ^ X i ” X 12 )cos a x + “ Z i2^ sin a ^ 

Z 1 “ (Z i " z i2 )cos a l " (X i " Z 12 )sin °1 


■ tan 



i 


1 , 2,. ...12 


(47) 

(48) 


For the transformation of points 12 through 23 to those in the X Y Z system, 
the following equations may be used: 


where 


x" 

X i 


( x i - X 12 )cos ct 2 


Z ’i 


( z i - Z 12 ) cos a 2 


a 2 - tan 


-1 


Z 12 ~ Z 23 


X 12 “ X 23 


+ (Z ± - Z 12 )sin a 2 
- (X ± - X 12 )sin a 2 

, i * 12,13,... ,23 


(49) 

(50) 


The (X^,Zp and (X^,Z^) points are curve-fitted and additional points are 
generated. Since these points must be transformed back to the reference 
system, equations (47) and (48) are solved for X^ and Z^ yielding 


Xj - Xj cos Oj - Zj sin + X J2 

Zj « Xj sin a 1 + Z* cos + Z 12 
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while equations (49) and (50) are solved for ^ and Z R yielding 

^ - XjJ cos a 2 - Zj[ sin a 2 + * 12 

Z k - Z£ sin o 2 + Z” cos a 2 + Z 12 

where the subscripts j and k denote appropriate, but different, point counters. 

Often during the preparation of data, errors are made either during the 
process of keypunching the data or through the lack of a clear visual or 
mental conception of the body's shape, which may be disastro-s to the poten- 
tial-flow calculations. These errors in the input are difficult to recognize 
until the data is plotted. 

The present method contains three techniques (options) to effectually 
correct these mistakes in a simple fashion. The first option modifies the 
data by explicitly changing the coordinates of individual originally-inputted 
body points through additionally-supplied coordinate information. By the 
second option, a body point P is replaced either by the average of the two 
points 1 and 2 (Figure 24a) of the M-line intersecting point P, the average 
of the two points 1 and 2 (Figure 24b) of the N-line intersecting point P, or 
the average of the four points 1, 2, 3, and 4 (Figure 24c) of the P-intersect- 
ing M- and N-lines. The third option uses the same point-identification 
arrangement as that of the second option, but replaces a body point P by 
linear interpolation (s) . Utilizations of these options are shown in Figure 
25 for a body assumed to be significantly in error at one point. 

The present method incorporates two more geometry-related options to aid 
in the refinement of the grid network. These options are important since the 
panel geometries of the inlet and exhaust sites must be relatively small 
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Figure 24. Coordinate information for averaging techniques 

compared to the overall body geometry to allow for the assumption that the 
opened-body pressure coefficients (see previous section) to remain essentially 
equal to those of the closed body. The first of these options augments the 
number of M-lines and/or N-lines by equally distributing the user-specified 
number of additional M- and/or N-lines along the arc length (assumed to be the 
linear distance) of every two successive M-lines for additional M-lines and 
along the perpendicular distance of every two successive N-lines for additional 
N-lines. The points on the additional M- and N-lines (equal in number and 
order to those on the original M- and N-lines) are calculated by cubic- 
polynomial spline-fit interpolations. This scheme, to be called the "equal- 
line augmentation", should be recognized to have the ability to increase the 
number of quadrilaterals (panels) manyfold with a corresponding decrease in 
the individual panel areas (Figure 26)* Although potential-flow calculations 
become more and more accurate as more and more panels are used to describe a 
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Figure 26. Example of equal-line augmentation scheme 
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body, core »u.t b. t.kon not to exc.d tho .v.lLblo oonput.r ototMd «p.bll- 
ities. 

For the above reason and that only a grid network refinement in the 
neighborhoods of the inlet and exhaust sites is usually desired, a second 
option has been included. Although it uses the same M- and N-line distri- 
bution techniques and point calculations of the first option, the second 
option augments the total number of M-lines and/or N-lines only by generating 
additional lines between any two specific M-lines and/or N-lines chosen by the 
user. Without a drastic increase in the total number of body panels, this 
scheme, to be called the "user-specified line augmentation", has the ability 
to generate additional lines on specific regions of the body while increasing 
the number of panels in these regions with individual areas smaller than the 
original panels (Figure 27). If the user prefers, this option may be used to 

duplicate the results _of_the first option. 

For a given body, aspects of both line augmentations may be used simul- 
taneously if the preferences of one do not override those of the other. That 
is, sometimes the equal-line augmentation scheme may be desired for M-lines 
(or N-lines) , while the user-supplied line-augmentation scheme may be desired 
for N-lines (or M-lines) . The implementation of these preferences are easily 
accomplished through the present method’s computational logic. 

Containing the coordinates of the points describing the original body 
(with the modifications invoked by the first three geometry options) and/or 
the final modified-body, cards may be punched in two different forms. The 
first form is compatible as„ input to the NCSU BODY potential-flow program 14] 
and that of the present report, while the second form is compatible as input 
to the NCSU PLOT program [4] for plotting complete configurations. 
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Given a set of body data (either input or_final), the grid-refinement 
program contains a plot option that allows for the generation of the necessary 
instructions for automatic plotting of the body and can be used to .draw three- 
view and oblique orthographic projections, as well as perspective and stereo- 
scopic projections. This option is a valuable tool for checking for inputted 
errors and for displaying the body modifications. Although examples of this 
plotting routine (originally written and programmed by Craidon, Reference 13) 
have been presented throughout this report, explicit exam ples are p resented in 
Figures 28, 29, 30, and 31. 

A geometry package has been developed by Halsey and Hess [Reference 14] 
which semi-automat ically panels isolated components, such as wings, fuselages, 
etc., of complete aircraft configurations by using one of several element 
(quadrilateral) distribution algorithms. Among other options for repaneling 
components, the Halsey-Hess program is similar to that of the present work 
since both programs allow the input of sparse coordinate data, use independent 
cubic-spline curve fits for interpolations, and provide many schemes for the 
point distributions on M- and N-lines. However, significant differences do 
exist. The curve-fit method used by Halsey and Hess does not insure continuity 
of the second derivative and thus is not a true cubic-spline fit in the usual 
sense. They claim their method gave consistently superior results to those of 
a true cubic spline. Their point-distribution schemes are also significantly 
different but are logical and perform to a fair degree of accuracy. In addi- 
tion, their geometry package contains a feature to calculate curves of inter- 
section among components and consequently repanel the regions of the components 
at the intersections. This feature should be extremely useful for design of 
complete configurations and the study of interference effects among the 
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Figure 29. Example of a oblique orthographic projection 
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component.. Sine, only Isolated coi.pon.nt. .r. .n.lyt.d by th. ptM.nt. work, 
this option is not available# 

For a geometry program to be .fflctent, It moat be .imple to use, h.ve 
accurate but expeditious algorithm., and preaant_ results logically for lmme 
dlate evaluation. While both programs are simple to use and have rapid and 
accurate algorithms, th. present program contains Its own bullt-ln plotting 
routines that allow the user to visual^ Inspect the point modifications and 
line augmentations. This feature eliminates the wasteful intermediate punch- 
ing* of cards as Input to separate plot programs, and therefore the overall 
time (program execution plus real time) to analyse body data Is significantly 
reduced.. The specialisation of the program to the specific analysis capability 
described In previous sections also reduces the else of the geometry program 
and may provide more expeditious execution as well. 


*The recording of data to magnetic tape or disk is common practice to avoid 
card punching. 
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The aforementioned methods have been implemented into computer programs. 
In addition to the boundary layer, the nonuniform propeller slipstream, and 
the cooling interior-mass-flow simulations, a program performing the potential 
flow calculations has been named FLOWBODY, while a program performing the grid 
modifications has been named GRIDPLOT. The listings of both programs - along 
with their user's instructions, sample input, and sample output - are provided 
in Volume 2. 

With logic easily adaptable to other computer facilities, both programs 
have been written in single-precision FORTRAN language. Single-precision 
arithmetic was chosen for two reasons. The first reason is that any smooth 
body represented by, an arrangement of finite-sized plane quadrilaterals is 


inevitably crude, while the second is that the increase in computational times 
(and therefore costs) for additional precision is undesirable. 
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DISCUSSION OF RESULTS 


An investigation has been made of means to reduce the drag of light air- 
craft fuselages and nacelles by recontouring the body and altering locations 
for the air intake and exhaust sites. Six bodies* were investigated with re- 
contouring performed on only one because of time limitations. A general pattern 
has emerged that satisfies intuition as well as the physical phenomena. As an 
aid in the discussion of the forthcoming results of the present method, certain 
aspects of theoretical and physical flow behavior will be reviewed briefly. 

The total or over-all drag on a body placed in a stream of fluid consists 
of skin friction (equal to the integral of all shearing stresses taken over 
the surface of the body) and of form or pressure drag (integral of components 
of the normal forces). Since it is important to determine the origin and mag- 
nitude of the skin friction and pressure drag, boundary-layer theory provides 
a means to explain these quantities in a rational fashion and also offers in- 
sight into what shape a body must take to minimize the total drag. 

In the immediate neighborhood of a solid wall, boundary- layer flow (laminar 
or turbulent) may, under certain conditions, become reversed causing the exter- 
nal flow to separate from the wall. Accompanied by the formation of eddies 
and large energy losses in the wake of the body, this separation phenomenon 
changes the pressure distribution which differs markedly from that in a friction 
less flow stream. The pressures on the lea side of the body are those of the 

*Three -dimensional bodies with a X-Z plane of symmetry (plane parallel 
to the flow direction.) 
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wake and substantially balow thosa on tha windward aide of the body which exper- 
ience near-stagnation pressures. Integration of such a pressure distribution 
shows a large drag. Hence, a large pressure drag nay be explained by the exist- 
ence of a large pressure deviation. 

Schlichting [15] gives an excellent description of boundary-layer separa- 
tion by considering the flow about a blunt body, e.g., about a circular cylinder. 
In frictionless (perfect fluid) flow, the fluid particles are accelerated on 
the upstream half and decelerated on the downstream half. Hence, in accordance 
to Bernoulli's theorem, the pressure decreases along the upstream half, while it 
increases along the downstream half. So long as the boundary layer remains 
thin, the flow is frictionless and pressure is constantly being transformed 
into kinetic energy along the upstream half. On the downstream half, kinetic 
energy reverts back to pressure in proper amounts such that the pressure at the 
rearmost stagnation point exactly equals that of ambient conditions. Since the 
fluid particles in the boundary layer remain under the influence of the same 
pressure field as that prevailing externally, the pressure distributions are 
ideal and the body does not experience a pressure drag. 

When viscosity (friction) is introduced, the particles expend so much of 
their kinetic energy during their travel along the upstream half of the cylinder 
that the remaining kinetic energy is sufficient for the particles to overcome 
the adverse pressure gradient on the cylinder's downstream half. Under this 
condition, the fluid motion within the boundary layer in the vicinity of the 
wall is eventually stopped and reversed by the external pressure field, result- 
ing in separated flow. The wide separation of the streamlines behind the cyl- 
inder and their failure to merge smoothly indicate that the flow is not of the 
constant-energy type, and therefore Bernoulli's equation is unsuitable for the 
determination of pressures on the surface in the wake. Within the wake region 
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behind the cylinder, a preaaure-diatribution curve clearly depicts suction or 
lea s-than-amb lent pressure. Because of this pressure difference, a pressure 
drag occurs. 

The shape of the body also controls the nature of the boundary-layer flow, 
i.e,, whether the flow is laminar or turbulent. For minimum skin-friction 
drag, it is necessary or logical to maintain a laminar boundar over the greatest 
possible distance of the body. Although efforts have resulted in a reduction 
of skin-friction drag, attempts to maintain laminar flow along too great a sur- 
face distance can produce large increases in the over-all drag from high pres- 
sure drag. Since pressure drag can be many times greater than the skin-friction 
drag and flow separation should be delayed as long as possible, a turbulent 
boundary layer is desirable. Even though a turbulent layer increases the skin 
friction, a turbulent boundary layer is more stable than a laminar one and is 
capable of maintaining itself longer under an adverse pressure gradient because 
of its greater capacity to remove energy from the freestream. Consequently, 
the flow adheres to the surface for a greater distance and, in turn, delays 
separation. The increase in skin-friction drag is usually much less than the 
reduction in pressure drag. 

Since the danger of boundary- layer separation always exists in regions of 
adverse pressure gradients with its likelihood of occurrence increasing for 
bodies with sharp or steep pressure curves (for instance, bodies with blunt 
ends), attempts should be made to seek a more streamlined body shape. The 
experimental pressure distributions for streamlined bodies differ so little 
trom those predicted for frictionless flow that the pressure rise in the down- 
stream direction is sufficiently gradual so that there is virtually no separa- 
tion. Consequently, the pressure drag is small enough that the total drag 
consists mainly of skin-friction drag. 
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When the boundary- layer equationa are integrated, the velocity distribution 
can be obtained, and the location of the point of flow separation can be deter- 
mined. This, in turn, allows the calculation of the viscous or skin-friction 
drag around the surface of the body by a simple process (in theory) of integrat- 
ing the shearing stress imposed by the fluid onto the surface. Even in cases 
where there is no_separation, the skin friction must be augmented by pressure 
drag because the boundary layer displaces the external potential* flow by a dis- 
tance equal to the displacement thickness. The pressure distributions are there 
fore changed even in the absence of flow separation, and the resultant of these 
pressure forces in the flow direction is no longer zero. When applied to actual 
problems, this method and the present method** allow the calculation of both 
the skin-friction drag, and the pressure drag.. 

Most fuselages and nacelles may be categorized as blunt or bluff bodies. 

For .this reason, their drag is almost entirely due to pressure deviations (from 
ambient conditions) and, therefore, is generally greater than that of stream- 
lined bodies of similar shape. In addition to rough surface textures, actual 
fuselages and nacelles experience interference effects from appendages and pro- 
tuberances that, normally, adversely affect the drag. Since few appendages 
(wings, landing gear, etc.) can be removed, protuberances such as protruding 
engine exhaust pipes, projecting rivet heads, improperly-designed cowl. flaps 
and canopies, etc. can be modified to provide the cleanest possible design 
without being impractical from the operational, maintenance, and financial 
points of view. 

*Flow external to the boundary layer is considered to be inviscid and 
therefore potential. 

**The reader should review the present report’s section on boundary-layer 
simulation in an inviscid- flow field. 
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The obvious complexity of actual flow fields about real aircraft bodies 
makes it necessary to investigate the influence of nomnifom* (in contrast to 
uniform) flows upon the over~all drag computations of the present methods 
Following the recommendations of a previous section for a crude approximation 
of a propeller slipstream, systems of ring vortices were placed about bodies in 
a systematic fashion so that there diameters, strengths, and locations reflected 
the presence of the bodies and a specified power into the airstream> These 
vortex systems induced velocities from approximately five percent of the free- 
stream velocity onto the body near the nose to about ten percent near the tail 
for a specified power of 100.66 kilowatts. 

The six bodies were tested with both uniform and nonuniform onset flow 
fields. The results are given in Figures 32, 33, 34, 35, 36, and 37 for the 
uniform flows, while those for the nonuniform flows are given in Figures 38, 

39, 40, .41,_42, and 43. As shown by the figures, five of the six bodies with 
a nonuniform flow field experienced reductions in skin-friction drag along with 
substantial increases in the pressure or form drag, resulting in larger total 
drags. The_other body also experienced an increase in pressure drag but was 
subdued by a substantial reduction in skin-friction drag to yield a lower total 
drag. A closer examination revealed that, possibly, the skin friction decreased 
because of the higher local Reynold's numbers at or near the body surfaces. 
Although the skin-friction drag was initially expected to Increase because of 
the additional "scrubbing" of the flow upon the bodies, a 30% increase in 
velocity may be seen to reduce the friction drag by 7% or more. Thus, these 
friction-drag reductions were acceptable. The pressure drag increased because 
of the higher velocities over the surfaces that significantly modified the 

*A nonuniform flow is defined as that onset flow with velocities differing 
from a constant value, whereas, an uniform flow is that onset flow with * con- 
stant velocity. 
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pressure distributions. The difference between the analytical method and 
physical phenonemon must also be noted. The analytical method automatically 
assumes flow separation near the rear of the body, while in actuality the flow 

may not separate because of the more energetic flow. 

The drag calculations of the ATLIT aircraft fuselage (Figure 4) and nacelle 
(Figure 42) were compared with the experimental drag findings by Holmes [l6] of 
a drag build-up for the complete aircraft. Holmes estimated the total aircraft 
drag coefficient to be 0.044. As a basis for comparison, the drag contributions 
[16] arising from the wing, the horizontal tail, and the vertical tail at zero 
angle of attack were accepted and used with the present method's prediction of 
those arising from the fuselage and nacelles to produce a total drag coefficient 
of 0.045. This prediction represents a significant improvement over previously- 
existing methods (including the 0.0358 value from earlier work of Smetana and 
Fox [17]) as well as portraying the present scheme as an useful tool for drag 

estimations. 

A "fat" nacelle (Figure 37 or 38) was contrived for the purpose of inves- 
tigating whether such a configuration might yield a lower drag for the ATLIT 
airplane. Initial calculations for the fat nacelle with an uniform onset flow 
field (Figure 37) showed a skin-friction drag coefficient of 0.00356 and pres- 
sure drag coefficient of 0.00317 to yield an over-all drag coefficient of 0.00673 
Initial calculations for the ATLIT nacelle with an uniform onset flow field 
(JEigure 36) showed a friction drag coefficient of 0.00308 and a pressure drag 
coefficient of 0.00414 to yield an over-all drag coefficient of 0,00722. Where- 
as, the friction drag for the fat nacelle was greater by 15.6%, a substantial 
decrease of 23.4% in the pressure drag was seen to cause a 6.8% reduction in 
the total drag. A nonuniform onset flow field was imposed onto the same two 
bodies to yield quite different results. For the fat nacelle (Figure 43), the 
skin-friction drag decreased by 9.2% to yield a skin-friction drag coefficient 
of 0.00323 and the pressure drag increased by 177% to yield a pressure drag 
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coefficient of 0.00879 for a total drag coefficient of 0.01203. For the ATLIT 
nacelle (Figure 42), the friction drag decreased by 7.5% to obtain a skin- 
friction drag coefficient of 0.00285 and. the pressure drag increased by 84% to 
obtain a pressure drag coefficient of 0.00762 for a total drag coefficient of 
0.01046. Obviously, the nonuniform flow had a significant effect since the 
total drag for the fat nacelle increased by 78.7% and that- for the ATLIT nacelle 
increased by 44.9%. Although the smaller friction drags were acceptable, the 
extremely large increases in the pressure drags were not. Since these nacelles 
were wider and flatter than most bodies considered, the concept and utilization 
of. the equivalent circular body (discussed in an earlier section) for the in- 
duced velocities by the propeller's slipstream onto the body surfaces were not 
suitable. For these nacelle cases and other possible cases of wide flat bodies, 
more accurate and realistic results would have been obtained with the use of 
an equivalent ellipsoidal body concept where the circular vortex rings are re- 
placed with elliptioally- shaped vortices. In addition to the unsuitability . 
of the circular body scheme, the present tests were conducted with a propeller 
diameter less than the maximum cross-sectional dimension of the nacelles. 
Consequently, the propeller forced the air flow onto (rather than about and. 
along ) the surfaces, causing an appreciable variation in the velocity and pres- 
sure distributions and, in turn, resulting in a large increase in pressure drag. 
Specification of larger propeller diameters should lessen the increase to more 
realistic values. 

Since all of the bodies were modeled by an orderly grouping of a large 
number of plane quadrilaterals (panels), it was therefore inevitable that such 
u modeling scheme produced only approximate representations to the actual 
bodies. For this reason, minor body recontouring was performed to study the 
severity of these crude models. The fuselage of a Cessna 182 airplane was 
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chosen for this test. The body was recontoured by the GRIDPLOT program, while 
the drag computations were effected by the FLOWBODY program. Both programs are 
contained in Volume 11 of this work. 

The study revealed that the total drag did not change appreciatively with 
typical calculations leading to drag reductions on the order of only one per- 
cent. In some cases, the total drag actually increased because of greater pres- 
sure drag attributed to larger adverse pressure gradients on the recontoured 
regions. In all cases the change in friction drag was insignificant. As an 
explicit- example, the upper rearward cabin section of a Cessna 182 fuselage 
(Figure 38) was slightly recontoured-. or "smoothed* 1 (Figure 39). Although dif- 
ficult to see these modifications because of the scale and orientation of the 
drawings, calculations yielded a 0.42% total drag reduction with approximately 
88% of this change arising from a decrease in the pressure drag. These results 
are encouraging since it is believed that further body modifications may reduce 
the drag even more. 

A different situation occurred when more severe shape changes were imposed. 
The nose of the representation of the Cessna fuselage was severed or flattened 
to better approximate the actual fuselage (Figure 40). When compared with the 
results for the original fuselage (Figure 38), this action produced an increase 
of 11.3% in the pressure drag and a slight reduction of 0.65% in the skin- 
friction drag. Like the case of a flat, plate nearly perpendicular to the flow 
direction, the increase in the pressure drag was entirely due to the disruption 
of the external flow field by the pressure of the blunt nose. In other words, 
the pressures on the blunt-nose panels of the fuselage approached the stagna- 
tion pressure, while the pressures on the "same" panels of the original fuselage 
were far less than the stagnation value. With the pressure distributions there- 
fore changed, integration of the pressures over the entire surface produced a 
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net increase in the pressure dreg. Hie skin* friction drag decreased primarily 
because of a decrease in the surface area. 

Although modifications to body shape may be seen to affect the drag more 
or less, the main purpose of this work is to determine suitable locations for 
the intake and exhaust sites of the air mass flow for effective cooling and 
minimum adverse disruption of the external flow field. Intuitively, the loca- 
tion and size of these sites should be dependent on or influenced by the power 
setting and the magnitude of the effective orifice area of the engine. 

Approximately 20 computer runs of the FL0WB0DY program were performed on 
the blunt-nose Cessna fuselage. For the fuselage .with the specification of the 
identical intake and exhaust sites and the same initial effective orifice area 
(Figures 44, 45, and 46), the total drag coefficient C D was found to vary 
linearly with power (Figure 47). The validity or pausibility of these drag 
coefficients may be seen by comparing the present results with those from wind- 
tunnel tests of a light single-engine aircraft [Reference 18] and a light twin- 
engine aircraft [Reference 19]. The power settings were converted to thrust 
coefficients, and the drag coefficients of the present method were plotted along 
with those from wind-tunnel tests corresponding to zero lift. It can be ob- 
served from Figure 48 that the present method yielded drag coefficients approx- 
imately 60-70% below those from the wind-tunnel tests. Whereas, the wind- 
tunnel tests considered the full aircraft configurations, the present drag- 
coefficient extractions pertained only to those of the fuselage and, therefore, 
should be only 30-40% as large. The curve representing the results of the 
single-engine aircraft's wind-tunnel tests is shown to behave remarkedly dif- 
ferent from that of the twin-engine aircraft and the present method's findings. 
Such behavior may be attributable to higher local Reynold's numbers or to small 
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or moderate flow separation. Further comparisons for similarly-shaped fuselages 
of References 9, 20, and 21 indicated that the extracted drag coefficients were.. 

indeed well within the range of typical values. 

Upon the designation of specific panels on the body as air intake ports or 
exhaust ports, the FLOWBODY program assumes that the pressure coefficients on 
the "opened" panels are equal to those existing on the same impermeable panels 
of the closed body. In actuality, the pressure coefficients may be quite dif- 
ferent if the respective panel areas are large in ^comparison to the total area 
of the body. Careful attention must be paid to satisfying this assumption if 
reasonable results are to be obtained. Since this analysis considers the 
interior of the body as essentially a "black box" containing an orifice (engine), 
a pressure drop normally occurs between the intake ports and the engine and, 
again, between the engine and the exhaust ports. Therefore, determination of 
two interior pressure coefficients is necessary to balance the mass flow pro- 
perly. Since che models fcr the program consist, of an arrangement of a large 
number of panels feach having its own pressure coefficient), the assumed pres- 
sure coefficients for the intake and exhaust porta are area-averages of those 
pressure. coefficients existing on the individual sites. 

As shown in Figure 49, an interesting and powerful result occurred when 
these pressure coefficients were plotted as functions of the effective orifice 
area, while the power and the intake and exhaust sites remained fixed. As the 
effective orifice area (EOA) was increased, the difference between the two 
interior pressure coefficients decreased. When EOA was increased to and beyond 
a value of 0.12 m 2 , the interior pressure coefficients had asymptotically 
approached a common value of 0.692. Therefore, a limiting value of EOA can be 
determined for each power setting and arrangement of intake and exhaust ports. 
Significant drag reductions were also observed as would be expected since the 
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engine is now effectively out of the flow peth. At these trees, cooling can be 
expected, to be minimal. As depicted in Figure 50 (extracted from Figures 51, 

52, 53, and 54), the total drag coefficient for the blunt-nose Cessna fuselage 
operating at the indicated test and freeatream conditions decreased to 0.01240 
(a 28.5% reduction in the drag of the same "closed" body) as the limiting EOA 
was reached. Since the_.mass flows entering and leaving the body are equal, 
sufficient air mass flow for effective engine cooling may be easily determined 
once the cooling requirements are known. It should also be noted that above 
the limiting EOA the drag no longer decreased but remained constant. Hence, a 
minimum-drag condition may be determined for each power setting and arrangement 
of the air intake and exhaust ports. 

Intuition would suggest that the limiting EOA should be related to the 
intake and exhaust areas. Since the intake and exhaust pressure coefficients 
are determined by area-averages which in turn limit the range of values for 
the interior pressure coefficients (see Figure 49), their magnitudes and signs 
wjuld indeed be functions of these areas. Upon this recognition, possible, 
"optimum" sites„.£or the air intake and exhaust may be chosen at the outset of 
the problem. 

Cowl air intake and exhaust sites must be suitably designed and located 
so that they provide adequate engine cooling and carburetor induction and pro- 
duce minimum pressure loss and favorable interactions with or minimum disruption 
to the external flows. Furthermore, they must provide good entrapment and 
extraction of the cooling air. The present method presumes that flows enter 
and exit the body normal to its surface and, consequently, must be of relatively 
low velocities to avoid severe disturbances to the external flow. It is advan- 
tageous to obtain prudent insight for judiciously selecting the proper sites. 
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A fundamental knowledge of fluid dynamica la effectually sufficient for this 
task. 

Several representative test cases were performed on the blunt-nose Cessna 
fuselage to yield the following general results: 

(a) Exhausting cooling air into concave regions of the fuselage at 
relatively low velocities results in significant decreases in 
over-all pressure drag, probably because it effectively increases 
body streamlining. Two typical regions where this is particularly 
effective are immediately upstream of the canopy and immediately 
downstream of the canopy. 

(b) The volumetric flow rate of cooling air passing through the fuse- 
lage should not exceed the amount required for cooling because 
significant drag increments are associated with excesses. Since 
the flow rate depends upon the flight conditions, provisions 
should be available for controlling this flow rate. This sug- 
gests the use of flush, sliding cowl flaps which efficiently 

r agulates the flow by varying the area of the exhaust port. 

(c) For good entrapment and extraction of the cooling air and 
minimum pressure loss, intakes should be located over regions 
with pressures near that of stagnation and exhausts should be 
located in low-pressure concave regions such that the exhaust 
flow produces a minimum disruption to or a favorable inter- 
action with the external flow. 

(d) Sharp edges, bluff frontal surfaces, and short afterbodies 
should be avoided. These are found to produce significant 
regions of flow separation and therefore increased drags. 
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The results of the 20 test esses conducted on the blunt-nose Cessns 
fuselage were insufficient to establish the ’’best" minimum-drag condltion(s) 
for a particular power setting. The majority of the cases were performed to 
observe the behavior of the results, and thus only a small number remained 
for further verifications and comparisons. Although limited in scope, three 
additional fuselage shapes were investigated to yield results consistent with 
those of the blunt-nose Cessna. Without much fear of contradiction, the re- 
peated occurrence of consistent patterns was perceived as generalizations 
applicable for a wide range of light aircraft. 
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Body Description: Cessna 182 Fuselage Model 

Flow Description: Uniform 



Test Conditions: 

Freestreara velocity (VINF), m/sec ••••••••• 48.768 

Freestreara density (ROE), kg/tn^ 1,2252 

Kinematic viscosity (VO), vcr/ sec 1.486 x 10 

Reference area (REFA), 16.165 

Developed engine power (DEP), kW 0.0 

Reynold's number 24,420,000 

Analysis Results: 

Pressure C T (w/o wake-body) .................. 0.00446 

L 

Pressure Cp (w/o wake-body) 0.00488 

Friction C Q 0.00780 

Pressure Cp (with wake-body) 0.00444 

Pressure Cp (with wake-body) 0,00491 


Total Body C L - 0.00444 
Total Body Cp ■ 0.01270 


Figure 32: Results of Cessna 182 fuselage model with uniform flow 

(power-off) . 
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Body Description: Recontoured Cessna 182 Fuselage Model 
Flow Descri ption: Unifora 



Test Conditions: 

Freest ream velocity (VINF), m/sec 48.768 

Freestream density (ROE), kg/m^ 1.2252 

Kinematic viscosity (VO), nr/ sec 1.486 ¥ 10 

Reference area (REFA), nr 16.165 

Developed engine power (DEP), kW 0.0 

Reynold's number ... 24,420,000 

Analysis Results: 

Pressure (w/o wake-body) C. 00425 

Pressure (w/o wake-body) 0.00507- 

Friction 0^ 0.00783 

Pressure (with wake-body) 0.00423 

Pressure C^ (with wake-body) 0.00509 


Total Body C^ ■ 0.00423 
Total Body C^ *» 0,01291 

Figure 33: Results of recontoured Cessna 182 fuselage model with 

uniform flow (power-off). 


83 


nntfilNAL PAGE IS 


Body Description: Blunt-nose Cessna 182 Fuselage Model 
Flow Description: Uniform 



Test Conditions: 

Freest ream velocity (VINF), m/sec 
Freestream density (ROE), kg/m • 
Kinematic viscosity (VO). m z /sec 
Reference area (REFA) , m ....... 

Developed engine power (DEP) , kW 
Reynold’s number - 

Analysis Results: 

Pressure (w/o wake-body) 
Pressure C D (w/o wake-body) 

Friction Cp 


Pressure C, (with wake-body) 

L 

Pressure Cp (with wake-body) 

Total Body C L - 0.00651 
Total Body C D - 0.01314 

Figure 34: Results of blunt-nose Cessna 182 fuselage model with 

uniform flow (power-off). 


0.00653 

0.00539 

0.00772 

0.00651 

0.00542 


48.768 
1.2252 
1.486 x 10' 
16.165 
0.0 

23,750,000 
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Body Description: ATLZT Fuselage Model 
Flow Description: Uniform 


\ 



Test Conditions: 


Freestreara velocity (VINF), m/sec 48.768 

Frees tream density (ROE), kg/m^ .*.!!!! 1.2252 

Kinematic viscosity (VO). mV sec ,* 1.486 x 10 

Reference area (REFA), m* 14.40 

Developed engine power (DEP), kW 0.0 

Reynold’s number 28,250,000 


Analysis Results: 

Pressure C L (w/o wake-body) 0.00329 

Pressure C D (w/o wake-body) 0.00222 

Friction C D 0.01057. 

Pressure C L (with wake-body) 0.00300 

Pressure Cp (with wake-body) 0.00242 


Total Body » 0.00300 

Total Body * 0.01299 


Figure 35: Results of ATLIT fuselage model with uniform flow (power-off), 
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Body Description: ATLIT Nacelle Model 
Flow Description: Uniform 



Test Conditions: 

Freest ream velocity (VINF), m/sec 
Freestream density (ROE), kg/m^ . 
Kinematic viscosity (VO). m z /sec 

Reference area (REFA), m z 

Developed engine power (DEP), kW 
Reynold’s number 


Analysis Results: 


Pressure C. (w/o waker-body) 0.00229 

L 

Pressure Cp (w/o wake-body) 0.00391 

Friction Cp 0,00308 

Pressure. (with wake-body) .. 0.00235 

Pressure Cp (with wake-body) 0,00414 


Total Body C L - 0.00235 
Total Body Cp ■ 0,00722 


48.768 
1.2252 
1.486 x 10 
, 14.40 
, 0.0 

, 9,667,000 


Figure 36: Results of ATLIT nacelle model with uniform flow (power-off). 


Test Conditions: 


Freest ream velocity (VINF), m/sec ..... 48.768 

Freestream density (ROE), kg/m^ 1.2252 

Kinematic viscosity (VO), m"/sec 1.486 x 10" 5 

Reference area (REFA), m* 14.40 

Developed engine power (DEP) , kW 0.0 

Reynold’s number 9,628,000 


Analysis Results: 

Pressure C L (w/o wake-body) 0.00228 

Pressure C D (w/o wake-body) 0.00331 

Friction C D 0.00356 

Pressure C L (with wake-body) 0.00240 

Pressure C D (with wake-body) 0.00317 


Total Body *• 0,00240 
Total Body Cp « 0,00673 


Figure 37: Results of a fat nacelle model with uniform flow (power-off) 




Body Description: Cessna 182 Fuselage Model 
Flow Description: Nonuniform 
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Test Conditions: 

Freest ream velocity (VINF), m/sec 
Freestream density (ROE), kg/m^ . 
Kinematic viscosity (VO), nr/sec 

Reference area (RE FA) , nr 

Developed engine power (DEP), kW 
Reynold's number 


48.768 
1.2252 
1.486 x 10" 
16.165 
100.66 
24,420,000 


Analysis Results: . 

Pressure (w/o. wake-body) 0.00358 

Pressure Cp (w/o wake-body) 0.00930 

Friction 0.00699 

Pressure (with wake-body) 0,00354 

Pressure Cp (with wake-body) 0,00934 


Total Body C L - 0.00354 
Total Body Cp ■ 0,01633 


Figure 38: Results of Cessna 182 fuselage model with nonuniform 

flow (power-on) . 
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Test Conditions: 


Freestream velocity (VINF), m/sec 
Freestream density (ROE), kg/m’ . 
Kinematic viscosity (VO), nr /sec 

Reference area (REFA), nr 

Developed engine power (DEP), kW 
Reynold's number 


48.768 
1.2252 
1.486 x 10 
16.165 
100.66 
24,420,000 


Analysis Results: 

Pressure C T (w/o wake-body) .. • ••• 0.00330 

L 

Pressure Cp (w/o wake-body) 0,00934 

Friction Cp 0.00699 

Pressure C L (with wake-body) 0.00326 

Pressure Cp (with wake-body) 0.00926 


Total Body ■ 0.00326 
Total Body Cp * 0.01623 


Figure 39: Results of recontoured Cessna 182 fuselage model with 

nonuniform flow (power-on). 
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Body Description: Blunt-nose Cessna 182 Fuselage Model 
Flow Description: Nonuniform 



Test Conditions: 

Freest ream velocity (VINF), m/sec 48.768 

Frees tream density (ROE), kg/m’ 1,2252 

Kinematic viscosity (VO). m z /sec 1,486 x 10 

Reference area (REF*,), m 16,165 

Developed engine power (DEP) , kW ..., 100,66 

Reynold’s number 23,750,000 

Analysis Results: 

Pressure C T (w/o wake-body) 0.00514 

Li 

Pressure Cp (w/o wake-body) 0.01035 

Friction Cp 0.00694 

Pressure C t (with wake-body) 0.00510 

la 

Pressure Cp (with wake-body) 0.01040 


Total Body Cp » 0.00510 
Total Body Cp » 0.01734 


Figure 40: Results of blunt-nose Cessna 182 fuselage model with 

nonuniform flow (power-on). 
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Test Conditions: 


Freestream velocity (VINF), m/sec 48.768 

Freestream density (ROE), kg/m^ 1.2252 

Kinematic viscosity (VO), nr/sec 1.486 x 10 

Reference area (REFA), nr 14,40 

Developed engine power (DEP), kW 100.66 

Reynold's number 28,250,000 

Analysis Results:. 

Pressure (w/o wake-body) 0.00246 

Pressure (w/o wake-body) 0.00251 

Friction C D 0,00927 

Pressure (with wake-body) 0.00223 

Pressure (with wake-body) 0.00271 

Total Body ■ 0,00223 

Total Body ■ 0,01198 


Figure 41: Results of ATL1T fusealge model with nonuniform flow 

(power-on) . 
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Body Description: ATLIT Nacelle Model 
Flow Description: Nonuniform 



Test Conditions: 

Freest ream velocity (VINF), m/sec 
Freest ream density (ROE), kg/m^ . 
Kinematic viscosity (VO). m z /sec 
Reference area (REFA), m z ....... 

Developed engine power (DEP), kW 
Reynold’s number 


48.768 
1.2252 _ 5 

1.486 x 10 
14.40 
100.66 
9,667,000 


Analysis Results : 

Pressure C L (w/o wake-body) 

Pressure C D (w/o wake-body) 

Friction 

Pressure (with wake-body) . . . . . 
Pressure C D (with wake-body) 

Total Body C L ■ 0.00189 
Total Body Cp - 0.01046 


0.00182 

0.00745 

0.00285 

0.00189 

0.00762 


Figure 42: Results of ATLIT nacelle model with nonuniform flow 

(power-on) . 
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Test Conditions: 


Freestreara velocity (VINF), m/sec 
Freest ream density (ROE), ke/m • 
Kinematic viacosity (VO), nr/sec 
Reference area (REFA), m 
Developed engine power (DEP), kW 
Reynold’s number 


48.768 
1.2252 
1.486. x 10 
14.40 
100.66 
9,628,000 


Analysis Results: 


Pressure 



0.00187 

C L 

yw/o waKe—Doay ) #•.•••••••« 

0.00894 

Pressure 

C D 

p 

yw/o waKC—ooay^ ••••*•***' 

0.00323 


Friction 

Pressure 

C D 

p 


.... 0.00199 

C L 

ywitn waK-e-ooay/ •••*••*• 

0.00879 

Pressure 

C D 

ywitn waKr"ouoyy *••••••• 



Total Body C L - 0.00199 
Total Body Cp » 0,01203 


Figure 43: Results of fat nacelle model with nonuniform flow 

(power-on) . 
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Body Description: Blunt-nose Cessna 182 Fuselage Model ... 

Flow Description: Nonuniform 



(Shaded areas denote ports for air intake and exhaust) 


Test Conditions: 


Freestream velocity (VINF) , m/sec 

Freestream density (ROE), kg/m 

Freestream temperature (TINF) , °K 

Kinematic viscosity (VO), m 2 /sec 

Constant-pressure specific heat of air (CPHA) , kJ/kg K 

Reference area (REFA) , m 2 

Heat of combustion of fuel (HVF) , kJ/kg 

Specific fuel consumption (SFC) , kg/kJ 

Developed engine power (DEP) , kW 

Initial effective orifice area (EOA) , m 2 

Reynold's number 


48.768 
1.2252 
288.86 
1.486 X 10' 
0.3095 
16.165 
41787.54 
1.5561 X 10 
74.563 
0.04645 
23,750,000 


Specification of 

non- 

zero 

normal velocities 


Input Code: MSI 

MS 2 

NSl 

NS2 

Panel No. 

Type 

1 

2 

1 

2 

1 

Intake 

3 

4 

1 

2 

5 

Intake 

4 

5 

1 

2 

7 

Intake 

5 

6 

1 

2 

9 

Intake 

6 

7 

1 

2 

11 

I (take 

7 

8 

1 

2 

13 

Intake 

8 

9 

1 

2 

15 

Intake 

16 

17 

5 

6 

111 

Exhaust 

17 

18 

5 

6 

113 

Exhaust 

18 

19 

5 

6 

115 

Exhaust 

19 

20 

5 

. 6 

117 

Exhaust 

20 

21 

5 

6 

119 

Exhaust 


Figure 44: Results of blunt-nose Cessna 182 fuselage model (nonuniform 

flow, 74.563 kW power, intake and exhaust ports). 


ORIGINAL PAGEMg 

OF POOR QUALITY 


Analysis Results: 




(Average C p ) lnlQt ■ 0.72902 

(AVSC.S8 C p ) oxhauat 

• ~0- 

2 

Inlet area, m 

- 0.07158 

Exhaust area, m 2 

• -0 

V* x 

- 0.681 

Final E0A, m 2 

• 0 


c - 0.64605 




p x 




c . 0.35667 




P 

r xx 



Pressure 

(w/o wake-body) 



Pressure 

Cp (w/o wake-body) 

0,00656 


Friction 

C D 

0.00666 


Pressure 

C. (with wake-body) . . . . 
L 

0.00273 


Pressure 

C p (with wake-body) .... 

0.00660 



Total Body C L - 0.00273 
Total Body C D - 0.01326 


Figure 44: Concluded 
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Body Description $ -Blunt-nose Cesanal82-Fuselage Model 
Flow Description: Nonunifom 



(Shaded areas denote ports for air intake and exhaust) 
Test Conditions : 


Freest ream velocity (VINF) , m/sec 48.768 

Freestream density (ROE), kg/m^ 1.2252 

Freestream temperature (TINF), °K 288.86 

Kinematic viscosity (VO), m^/sec o’**** * 10- 

Constant-pressure specific heat of air (CPHA) , kJ/kg K ... 0.3095 

Reference area (REFA) , 16.165 

Heat of combustion of fuel (HVF) , kJ/kg 41787.54 

Specific fuel consumption (SFC) , kg/kJ 1.5561 X 10 

Developed engine power (DEP) , kW 100.66 

Initial effective orifice area (EOA) , m2 0.04645 

Reynold's number 23,750,000 


Specification of 

non- 

zero 

normal 

velocities 


Input Code : MSI 

MS2 

NSl 

NS2 

Panel No. 

Type 

1 

2 

1 

2 

1 

Intake 

3 

4 

1 

2 

5 

Intake 

4 

5 

1 

2 

7 

Intake 

5 

6 

1 

2 

9 

Intake 

6 

7 

1 

2 

11 

Intake 

7 

8 

1 

2 

13 

Intake 

8 

9 

1 

2 

15 

Intake 

16 

17 

5 

6 

111. 

Exhaust 

17 

18 

5 

6 

113 

Exhaust 

18 

19 

5 

6 

115 

Exhaust 

19 

20 

5 

6 

117 

Exhaust 

20 

21 

5 

6 

119 

Exhaust 
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Figure 45: Results of blunt-nose Cessna fuselage model (nonunifom 

flow, 100,66 kW power, Intake and exhaust ports). 
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Body Doacription: Blunt-nose Cessnal82 Fuselage Model 
Flow Description:. Nonuniform 



(Shaded areas denote ports for air intake and exhaust) 
Test Conditions: 


Freest ream velocity (VINF) , m/sec 48.768 

Freestream density (ROE), kg/m^ 1.2252 

Freestream temperature (TINF), °K 288.86 

Kinematic viscosity (VO), m^/sec 1.486 X 10”^ 

Constant-pressure specific heat of air (CPHA) , kJ/kg°K ... 0.3095 

Reference area (REFA), m2 16.165 

Heat of combustion of fuel (HVF) , kJ/kg 41787.54 

Specific fuel consumption (SFC), kg/kJ 1.5561 X 10"*' 

Developed engine power (DEP) , kW 223.69 

Initial effective orifice area (EOA) , m2 0.04645 

Reynold's number 23,750,000 


Specification of 

non- 

zero 

normal velocities 


Input Code: MSI 

MS2 

MSI 

NS2 

Panel No. 

Type 

1 

2 

1 

2 

1 

Intake 

3 

4 

1 

2 

5 

Intake 

4 

5 

1 

2 

7 

Intake 

5 

6 

1 

2 

9 

Intake 

6 

7 

1 

2 

11 

Intake 

7 

8 

1 

2 

13 

Intake 

8 

9 

1 

2 

15 

Intake 

16 

17 

5 

6 

111 

Exhaust 

17 

18 

5 

6 

113 

Exhaust 

18 

19 

5 

6 

115 

Exhaust 

19 

20 

5 

6 

117 

Exhaust 

20 

21 

5 

6 

119 

Exhaust 


Figure 46: Results of blunt-nose Cessna 182 fuselage model (nonuniform 

flow, 223.69 kW power, intake and exhaust ports). 
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ORIGINAL PAGE IS 
OF POOR QUALITY 

Analysts Results: 

(Average C p ) inlet ■ 0,6(1655 -(Average 0 p ) exhau8t _ 

Inlet area, m 2 - °* 07158 - Exhaust area, m 2 

T » 0.522 Final EOA, m 2 

r » 0.44380 

°p K 

r » -0.29653 

p 

r xx 

Pressure C T (w/o wake-body) 0,00051 

JL 

Pressure C D (w/o wakerbody) 0,02648 

Friction C D 0.00768 

Pressure C T (with wake-body) 0.00037 

L 

Pressure (with wake-body) 0,02658 




Total Body C L ® 0.00037 

Total Body C D a 0.03427 


Figure 46 : Concluded 
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Power, kW 

Figure 47: Variation of total drag coefficient Cd with power for blunt-nose 

Cessna fuselage model with fixed intake and exhaust sites. 


NASA TN D-6238 (Twin Engine) 


NASA TN D-5857 (Single Engine) 


N 


Present Method’s Predictions 


0 M !o4 " 0 6 .08 AO AZ ^4 AO .18 .20 

Ct 


Figure 48: Comparison of drag coefficients Cj) between prediction and 

wind-tunnel tests* 
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(Shaded areas denote ports for air intake and exhaust) 


Test Conditions: 


Freestream velocity (VINF) , m/sec . 

Freestream density (ROE), kg/m-^ 

Freestream temperature (TINF), °K 

Kinematic viscosity (VO), tn 2 /sec *•••*••**•**•* 

Constant-pressure specific heat of air (CPHA) , kJ/kg K . . . 

Reference area (REFA), m 2 

Heat of combustion of fuel (HVF) , kJ/kg 

Specific fuel consumption (SFC), kg/kJ 

Developed engine power (DEP) , kW 

Initial effective orifice area (EOA) , m 2 

Reynold's number 


48.768 

1.2252 

288.86 s 
1.486 X 10" 3 
0.3095 
16.165 

41787.54 . 

1.5561 X 10" 4 
74.563 
0.02322 
23,750,000 


Specification of 

non- 

zero 

normal velocities 


Input Code: MSI 

MS2 

NS1 

NS2 

Pauel No. 

Type 

1 

2 

1 

2 

1 

Intake 

3 

4 

1 

2 

5 

Intake 

4 

5 

1 

2 

7 

Intake 

5 

6 

1 

2 

9 

Intake 

6 

7 

1 

2 

11 

Intake 

7 

8 

1 

2 

13 

Intake 

8 

9 

1 

2 

15 

Intake 

16 

17 

5 

6 

111 

Exhaust 

17 

18 

5 

6 

113 

Exhaust 

18 

19 

5 

6 

115 

Exhaust 

19 

20 

5 

6 

117 

Exhaust 

20 

21 

5 

6 

119 

Exhaust 





Figure 51: Results of blunt-nose Cessna 182 fuselage model (nonuniform 

flow, 74.563 kW power, 0.023. m 2 initial EOA, intake and exhaust ports) 
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Analysis Results: 

(Average C p ) lnUt 
2 

Inlet area, m 

V*x 


■ 0.72902- 

- 0.07158 

- 0.542 
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(Average C p ) exhauBt 
Exhaust area, 

Final EOA, 



0.70938 

0.36477 


Pressure (w/o wake-body) 0.00384 

Pressure (w/o wake-body) 0.00787 


Friction C Q 


0.00665 


Pressure C T (with wake-body) 0.00381 

Pressure (with wake-body) 0.00791 


Total- Body C L ® 0.00381 

Total Body C D - 0.36477 




0.32620 

0.05109 

0.01709 


Figure 51: Concluded 
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Body Description: Blunt-nose Cessna. 182 Fuselage Model- 
Flow Description: Nonuniform 



(Shaded areas denote ports for air intake and exhaust) 
Test Conditions^ „ 


Freestream velocity (VINF) , m/sec 

Freestream density (ROE) , kg/m^ 

Frpestrcam tpmDpraf tirp (TTNF} , °K .......... .............. 

48.768 

1.2252 

odq ae. 

If i non. 1 ? t *! o vl Rons'! tv ( VO^ . tn^/sfir ......................... 

400 • OO - 

1.486 X lo" 3 

0.3095 

16.165 

41787.54 

1.5561 X 10" 

74.563 

0.04645 

23,750,000 

Constant-pressure specific heat of air (CPHA) , kJ/kg°K 

Rp f prpnnp arpn fREFA^ . m2 ................................ 

Heat of combustion of fuel (HVF), kJ/kg 

Specific fuel consumption (SFC) , kg/kJ 

Developed engine power (DEP) , kW 

Initial pffertivp orifice area ^EOA^ * m2 . ... 

Reynold's number 


Specification of 

non- 

zero 

normal velocities 


Input Code : MSI 

MS 2 

NS1 

NS2 

Panel No. 

Type 

1 

2 

1 

2 

1 

Intake 

3 

4 

1 

2 

5 

Intake 

4 

5 

1 

2 

7 

Intake 

5 

6 

1 

2 

9 

Intake 

6 

7 

1 

2 

11 

Intake 

7 

8 

1 

2 

13 

Intake 

8 

9 

1 

2. 

15 

Intake 

16 

17 

5 

6 

111 

Exhaust 

17 

18 

5 

6 

113 

Exhaust 

18 

19 

5 

6 

115 

Exhaust 

19 

20 

5 

6 

117 

Exhaust 

20 

21 

5 

6 

119 

Exhaust 
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Figure 52: Results of blunt-nose Cessna 182 fuselage model (nonuniform 

flow, 74,563 kW power, 0.046 ra z initial EOA, intake and exhaust ports) 
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Analysis Results: 

(Average y inlot 

- 0.72902 

- 

(Average C p ) exhaust 

» 0.32621 

2 

Inlet area, in 

- Oj, 07158 


Exhaust area, m 2 

* 0.05109 


» 0.643 


Final EOA, m 2 

- 0.03726 


n » 

o 

<** 

c p * 

r xx 

0.66847 

0.44504 




Pressure C L (w/o wake-body) . 0.00202 

Pressure (w/o wake-body) 0.00690 

Friction 0.00666 

Pressure (with wake-body) 0.00199 

Pressure. (with wake-body) 0.00694 


Total Body C t - 0.00199 

Total Body C D = 0.01360 


Figure 52: Concluded. 
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Body Description: Blunt-nose Cessna 182 Fuselage Model 

Flow Description: Nonuniform 



womrnmp 


(Shaded areas denote ports for air intake and exhaust) 


Test Conditions: . 


Freestream velocity (VINF) , m/ sec • • - * 

Freest ream density (ROE), kg/m-^. 

Freestream temperature (TINF), °K 

Kinematic viscosity (VO) , m 2 / sec • • • 

Constant-pressure specific heat of air (CPHA) > kJ/kg K • • • 

Reference area (REFA), m 2 

Heat of combustion of fuel (HVF), kJ/kg 

Specific fuel consumption (SFC) , kg/kJ * 

Developed engine power (DEP) , kW • • 

Initial effective orifice area (EOA), m 2 

Reynold's number 


48.768 
1.2252 
288.86 «. 
1.486 X lO -3 
0.3095 
16.165 
41787.54 
1.5561 x 10” 
74.563 
0.0929 
23,750,000 


Specification of 

non- 

zero 

normal velocities 

Type 

Input Code : MSI 

MS2 

NSl 

NS2 

Panel No. 

1 

2 

1 

2 

1 

Intake 

3 

4 

1 

2 

5 

Intake 

4 

5 

1 

2 

7 

Intake 

5 

6 

1 

2 

9 

Intake 

6 

7 

1 

2 

11 

Intake 

7 

8 

1 

2 

13 

Intake 

8 

9 

1 

2 

15 

Intake 

16 

17 

5 

6 

111 

Exhaust 

17 

18 

5 

6 

113 

Exhaust 

18 

19 

5 

6 

115 

Exhaust 

19 

20 

5 

6 

117 

Exhaust 

20 

21 

5 

6 

119 

Exhaust 
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Figure 53: Results of blunt-nose Cessna 182 fuselage model (nonuniform 

flow, 74.563 kW power, 0.093 m 2 initial EOA, intake and exhaust ports). 
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Analysis Results: 

(Average C ? ) inUt 
2 

Inlet area, m 

*Jh 


0.72902 

0.07158 

0.689 


(Average C p > exhaust 
Exhaust. .area, m 2 
Final EOA, m 2 


» 0.62370 
» 0.53290 


Pressure 


0*0037 

(w/o waK>e-D9ay; ■••••••• • 

0.00574 

Pressure C Q 

(w/o waKe-Doay; « 

0.00666 

Friction 



Pressure C 

li 


0.00033 

(with waK.e-Doay; 

0.00578 

Pressure 0^ 

(witn watve-Doay; ••••••• 



Total Body C L B 0.00033 
Total Body Cp = 0.01244 


Figure 53: Concluded. 
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Body Description: Blunt-nose Cessna 182 Fuselage Model 

Flow Description: Nonuniform 



(Shaded areas denote ports for air intake and exhaust) 
Test Conditions: 


Freestrearo velocity (VINF) , m/sec 

Freestream density (ROE), kg/m^ 

Freestream temperature (TINF), °K 

Kinematic viscosity (VO), m^/sec 

Constant-pressure specific heat of air (CPHA) , kJ/kg°K ... 

Reference area (REFA) , m 2 

Heat of combustion of fuel (HVF), kJ/kg 

Specific fuel consumption (SFC), kg/kJ 

Developed engine power (DEP) , kW 

Initial effective orifice area (EOA) , m2 

Reynold ' s number • • 


48.768 
1.2252 
288.86 . 
1.486 X 10" 3 
0.3095 
16.165 
41787.54 
1.5561 X 10“ 
74.563 
0.9290 
23,750,000 


Specification of 

non- 

zero 

normal velocities 


Input Code : MSI 

MS2 

MSI 

NS2 

Panel No. 

Type 

1 

2 

1 

2 

1 

Intake 

3 

4 

1 

2 

5 

Intake 

4 

5 

1 

2 

7 

Intake 

5 

6 

1 

2 

9 

Intake 

6 

7 

1 

2 

11 

Intake 

7 

8 

1 

2 

13 

Intake 

8 

9 

1 

2 

15 

Intake 

16 

17 

5 

6 

111 

Exhaust 

17 

18 

5 

6 

113 

Exhaust 

18 

19 

5 

6 

115 

Exhaust 

19 

20 

5 

6 

117 

Exhaust 

20 

21 

5 

6 

119 

Exhaust 




Figure 54: Results of blunt-nose Cessna 182 fuselage model (nonuniform 

flow, 74.563 kW power, 0.929 m^ initial EOA, intake and exhaust ports). 
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Analysis Results; 




(Averogo Vimec 

- 0.72902 

^Average C p > exhau8t _- 

0.32621 .. 

2 

Inlet area, in 

- 0.07158 

Exhaust area, m 2 ■ 

3,05109. 

V T i 

- 0.708 

c = 0.59344 

\ 

c « 0.59230 

p 

r xx 

Final EOA, m 2 

0.78153 


i 


\ ", 


1 


Pressure C T (w/o wake-body) 

L 

-0.00069 

Pressure C Q (w/o wake-body) 

0.00493 


0.00666 

r ric exon Vjj 


Pressure C T (with wake-body) ........ 

L 

-0.00073 

Pressure (with wake-body) 

0.00574 


Total Body C L ® -0.00073 
Total Body C Q - 0.01240 

Figure 54: Concluded. 
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CONCLUSIONS 


>■ * — - - "TrrrriTrrrci.. 

““tt:;: ... ...«•■ - — 

. j cizine of the cooling air 

28 . 5 % by body recontouring and proper placemen 

intakes and equate. ^ ^ ^ mII alrp l a „e using this 

3. The predicted total drag co _ 

1 2 3% large r than that obtained by Holmes [l6 us ng o 

approach was on y • . fuse la g e and nacelle drag pre- 

u • „«, Since the present work pertained to fuselag 

techniques, since uie y -irientical to 

dictions, drag contributed from other aircraft components uere made 

those of Holmes. _ „nrtices appears to 

1 1 ci-r flow bv systems ot ring 

4. Simulating the propeller flou by , 

st isolated bodies. The equivalent circula 
be adequate £or use with most an elHpsoidal- 

concept may be unsuitable Tor uide Hat bodies. In euch cases, 
body concept would be preferable. 
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